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I. Introduction 


The stopping power of a material is defined as the space-rate of 
loss of energy suffered by a particle of nuclear or atomic dimensions 
in its passage through the material. Since the time of Thomson and 
Rutherford the subject has been investigated both theoretically end 
experimentally in its many complexities, and it early provided a 
source of tests of atomic structure and mechanics, In recent years 
the interest of theoretical nuclear physicists in the various facets 
of the subject has waned, since it is largely a property of the extra- 
nuclear part of the atom. Nevertheless the subject is still a vital 
one in the experimental and applied side of nuclear physics. The 
need for good stopping power formulas and/or data for all types of 
nucleer particles, over the whole range of Moda and through all 
types of matter is almost too obvious to mention - 6.g., in cosmic 
ray studies, range determinations, shielding requirements, and so on, 

Unfortunately for those who need this information, stopping 
power formulas provided by theoreticians are in most cases unable to 
give accurate predictions without the use of semi-empirical constants 
adjusted to conform to previous experiment. This is due to the 
physical and mathematical complexities of the subject. By this date 
theory and experiment are sufficiently advanced thet for heavy par- 
ticles, such as protons or alpha particles, the stopping power can 
be predicted with accuracy for мй types of mtter usually encountered, 
provided the energy of the particles exceeds a few million electron- 
Tdv. However if one passes to the energy region lower than thís, the 


predictions became less accurate and even empirical adjustments are 
-/ 


έ 





difficult. In such regions experimental data must largely be relied 
on at present. Such data are still far from complete, even for 
stopping by basic elements. 

For direct measurement of stopping power within this energy 
region the Van de Graaff generator is well suited, because Sn 
range of accelerating potentials evailable and because of its ability 
to produce accelerated particles of precisely defined and accurately 
measured energies. For this dissertation the writer has elected to 
report on the experimental results of stopping power of various ele- 
ments for protons from the Van de Graaff generator at the Ohio State 
University. Specifically, measurements are herein reported of the 
stopping power for protons of the metals, copper and nickel, and of 
the gases, nitrogen, neon, argon, krypton, and xenon. The energy 
range over which these measurements are taken is from 400 to 1200 kev, 

Brief mention is made of the associated problem of "straggling" 
in stopping power, or the mean square deviation in energy loss per 
unit thickness of matter for individual protons. This phenomenon is 
a result of the fact that slowing dovm of particles is accomplished 
through random collisions with atoms of the stopping medium and is 
therefore subject to statistical variations. This factor is less 
important than that of average stopping power, and both theoretically 
and experimentally it has been studied less fully. Accurate analysis 
of straggling is just as difficult in theory as analysis of stopping 
power itself, and experimentally it is even more difficult, requiring 
generally a precision of measurement better than stopping power 


measurement by an order of magnitude. The techniaue employed by the 
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writer for determination of stopping does make it possible to deter- 
mine variations in same, under favorable experimental conditions. It 
Ф. been found possible to report stopping power stragsling for the 
metallic elements, copper and nickel, although precise accuracy is 
lacking. Such results are considered as a "by-product" of the major 


investigation into stopping power itself. 


A. Definitions and Units. 

Stopping power has been defined above; and unless otherwise 
specified or implied it refers to the average energy loss of a large 
number of similar particles of the same energy, per unit length of . 
path through the material traversed. It is dependent upon the tyne 
of particle, the energy (or velocity) of the particle, and the stop- 
ping medium. The concept of stopping power involves a continuous 
change in particle energy, and yet a specified value of stopping 
power is valid only for a fixed value of particle energy. ‘this in- 
dicates that stopping power must be recognized as a derivative, 

- dE/àx, itself a function of energy. 

If the stopving medium is composed of a single element, one can 
use the concept of atomic stopping power, defined as stopping power 
divided by the number of atoms ser unit cube of the stopping material, 
(This is sometimes called the "stopping cross-section" of the atom.) 
It may be considered as the average contribution of each atom of the 
elemental material toward slowing the incident particles. In symbols 
it is designsted as O” = = -&) e 


Theoretical formulas have been derived in the Gaussian (c.g.s.) 


= 





a 


system of units, so that stopping power is given in ergs/em. For 
experimental work - both in determination and use of stopring power - 
energy changes are more usually expressed in electron-volts or 
similear units, while thicknesses are conveniently expressed in mg/em", 
Thus stopping power data from experiment is usually reported in units 
of kev-cm” /me. The conversion factor between the two system of units 
is easily seen to be: 

Stopping power (ergs/cm) = 1.602 х 10-6 x Stopping power 

(kev-em"/mg) x P » (I-A~1) 
where 

2 = density of medium (g/emŠ) , 

For atomic stopping power the proper c.g.s. units sre erg-cm” 2 
Experimental conventions make preferable the units ev-em”, The 
conversion factor between the two systems is obviously: 

At, stopping power (erg-em°) a 1,602 x 10-15 x At. stopping 

power (ota? ) 3 (I-A-2) 

The conversion factor from OPRIMA power to atomic stopping 

power in the c.g.s. system is: 


(Stopping power) ¿pg [N 


(At. stopping pover) 
CES 


A 
23 
6.023 x 10 xo 


x (Stopping 





over I-A-3 
p — , ( ) 
where N is the number of atoms per cubic centimeter of stopping 
material, and A is the atomic weight of the stopping element, 
From the above equations it is easily deduced that in experimen- 


tal units the conversion factor from stonvping power to atomic stonning 


— ad. 
Ф 





power is: 


(At. stopning power) = A x (Stopping 
e 


SP 6.023 x 1017 
PR (1-Λ-4) 
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The concept of stopping power "straggling" is understood most 
easily by visualizing a statistical distribution of relative fre- 
quencies of stonping nower losses for individual particles, under 
identical conditions. If the standard deviation about the mean of 
this distribution is designated as SL, then f) becomes the mean 
square fluctuation , or the variance; and one best describes the 
straggling by = (NM) ¿e (In finite increments we may consider 
straggling as 0/ax , which is used rather than fL/Ax since, 
as we shall see later, the former is independent of thickness of 
medium while the latter is not.) This quantity may also be written, 
as shown by statistical theory (Hoel, 1947), as = ΚΕ - ay] ; 
where E in this expression indicates the energy of sn individual 
particle among a large population of sinilar nerticles incident on 
the stopping medium with the same energy. 

The units for E aa are evidently ergs” /em in the 4.@.8. 
system. In experimental work the units emploved ere kev^ -em* /ng e 
The conversion relation between the values of straggling as expressed 


in these two systems is: - 


e 
Е | N | cgs 7 101° = V^ А dx ( N ) exp. ° (I-A-5) 


B. Review of Stopping Power Theory. 


A complete review of stopping power theory would be almost im- 
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possible and certainly impractical for a report of this nature. It 

is necessary therefore to present only a brief outline of the theory. 
Much discussion can be eliminated by concentration only on those 
aspects bearing on our particular problem. For example, relativity 
considerations may be ignored, since we are interested only in a 
range of velocities well below that of light (v/c ~ 1/30). 
Furthermore, we need consider only the stopping of heavy particles 
such as protons, and primarily by the medium or heavy elements, 

In a broad sense there are two approaches to the problem, and 
each will be summarized in turn.* The first, exemplified by the 
works of N. Bohr and his followers, is classical or semi-clessical 
in spirit. This is advantageous for ease of understanding, mathe- 
matical simplicity, and flexibility in application to a wide variety 
of conditions. The most rigorous formulation of the subject has been 
achieved, however, by modern wave-mechanical methods, with their 
mathenatical complexities. Such methods may be difficult to ap- 
preciate intuitively, but they yield more accurate results in those 
realms where the besic assumptions are valid. 

In presenting the simpler theory we shall follow primarily the 
excellent treatment given by N. Bohr (1948). 

1. Classical Fundamentals. 
Let us consider the elastic collision between an incident 


particle (number one) and a second particle which is at rest in the 


1 
For the sake of clarity and brevity we are not following strictly a 


historical, or chronological, scheme of presentation. All major 
contributors to the theory however will be duly acknowledged where 
their ideas best fit into this summary. 
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laboratory frame of reference. Figure 1 shows the geometric rela- 
tionships existing between the velocities in the center-of-mass . 
system. 
E Y 
vo! ее 
ΤΗ 
Y-Vc 8 en -Va 


C eMe 


Figure T. 


In this figure, 


ql 


= velocity of particle l, laboratory system, before collision; 


«| 
" 


vel. of center-of-mass, laboratory system, before collision; 
x = vel. of particle 2, center-of-mass system, before collision; 
У-Ү = vel. of particle 1, center-of-mass system, before collision; 
6 » angle of deviation of each particle in the c.-m. system. 
(For direct collision, the particles sre always oppositely directed in 
the center-of-mass system, so that each particle undergoes the sane 
angular deviation in this systen.) 
Each particle behaves in the center-of-mass system as if it 
rebounds elastically from the immovable center of mass without loss of 


energy. Then, 


L^ = vel, of particle 2, center-of-mass system, after collision; 
and ^ а 

UM = |, aen . (I-B-1) 
Also, 





Y, = vel. of particle 2, laboratory system, after collision, 
= Yo + Y, » (I-B-2) 
We see that Vo forns the base of an isosceles triangle with 
vertex angle 8 . Therefore 
vo, = 2v, sin (0/2) . (1-B-3) 


From the definition of center-of-mass it is easily shown that 


m 
we —., , (I-B-4) 
Dy + Me 
where my and m, are masses of the two particles. Thus 
Yo = 2v—_ 1 sín (0/2) . (1-B-5) 
ἪΝ 


Now in the collision particle 2 is given energy at the expense 
of particle 1 (laboratory coordinates). Let T be this energy trans- 
ferred. 

T = m Vo" 9 (1-8-6) 
and from equation (I-B-5), 


2 
- τα 5: ү sin“ (0/2) 


(m +m)? 
(I-B-7) 


2 
2% 2 sin2(0/2) , 
mo 
where m, is the "reduced mass" and equals mym,/(m) + m,) » 
The maximm transfer of energy occurs on direct collision, when 
Ө equals 1800 , Let T, Pe this amount. Then 


2 
Ты = Zu y 


m = (I-B-8) 





and 
T = T, 8112(9/2) ° (I-B-9) 


If the law of attraction or repulsion between the two perticlés 
-8 





mE 





is Coulombian, one may derive the well known Rutherford scattering 


formula: , 
о 
do = [252 \ . све (0/2) - ἁω , (1-B-10) 
2 Mo ve 
where 


do = differential probability cross-section that the particle 
is scattered into solid angle, dw , measured in the 
center-of-mass coordinate system. 

dw = total solid angle corresponding to scattering directions 
between O and O+ de, 

з 217510 Ө de = 4Tsin(0/2) - cos(0/2) - 40 . 
(I-B-11) 
өр = charge of the first particle. 


9 charge of the second particle, 


Ú 
Substitution of the expression (I-B-11) into equation (I-B-10) gives: 


e) e> | 
åo- = 4T (23) - евв?(Ө/2) :eos(0/2) - a0 , (I-B-12) 
2m, v^ 


From equation (I-B-7) we obtain by differentiation: 
dT = T cos(0/2).csc(0o/2).d0 . (I-B-13) 


Then from the last two equations we obtain: 


2 2 
do s Ba 2) a 


mo v? тё 
(I-B-14) 
= p Sr 
me f 


where D is the quantity within the parentheses. 

We may use this result to obtain stopping power relations in a 
highly idealized case. Let us suppose that an incident particle 
undergoes a succession of energy-diminishing collisions in traversing 


a thin layer of matter in a direction normal to the layer, 
`< 





Let 
Ax 


N°’ 


thickness of layer of matter; 


mimber of particles per unit cube of the stopping 
medium that absorbs the energy. 

It is necessary to assume here that: 

(a) The particles of the stopping medium are free particles at rest, 
interacting with the incident part. through the Coulomb law of farce 
between charged particles. 

(b) The energy transfer in an individual collision is small compared 
with the total energy of the incident particle. This is true provided 
Te << 2 m v y Which according to equation (I-B-8) is true when 
m¿2/m Mm < 2. 

(ο) The medium layer is very thin so that v is essentially constant 
throughout the progress of the incident particle within the layer. 

Assumption (b) indicates that momentum change of the incident 
particle is small in a single collision, which permits one to state 
that the particle velocity is almost constant in direction as well as 
in magnitude. 

Now consider that we have made a series of hypothetical experi- 
ments in which we were able to follow each particle in an incident 
beam with common initial characteristics, and further that we could 
obtain and tabulate the data for each collision each particle under- 
goes with particles of the stopping material. Divide the possible 
values of energy loss in a single collision into equal intervals, 
with T corresponding to the average energy for the ‚th interval 


1 


and 04 being the probability that each collision will provide an 


- “2 








energy loss within the interval characterized by Ty: Since the 
stopping process is a statistical process because of the randomness 
of location of particles of the medium in relation to the path of 
some incident particle, we cannot always get in different trials the 
same number of collisions characterized by T4 . Define 
n, = number of collisions characterized by energy loss T4 , 
for the case of any one incident particle; 
πι = average number of collisions characterized by Ty , for 
the total population of incident particles in the beam, 
From the given definitions it is evident that 
w = NAI- . (1-В-15) 
In a practical case the number of collisions is large and 
individual "single shot probability" is small. In such catio the 


probability distribution function for n, can be closely approximated 


i 
by a Poisson distribution formula (Hoel, 1947): 


n = : 
Ey MAT, (I-B-16) 


having mean value w, as required and having standard deviation eaual 


to үтү Р 


Now if AE is the change in emergy of a particular particle 


i 


when passing through the given layer of the medium, 


-ΔΕ = δὴ . I-B-17 
A =, ¿Y ( ) 

The average value of energy loss 18 
-AE = = тҮ Wi ; (I-B-18) 


and the mean square fluctuation is by definition 


Ë E A 
f. = (Ag XE? 





= = тү? (n, - ye 


i 1 
= == 72 (Vw y? Kerns... (I-B-19) 
n i i 1 f 
If we go over to integrals and use equation (I-B-15) we obtain: 
- Тах 
У = NrÁx: T do- | (I-B-20) 
Tiin 
and 
2 Tnax 
(1 = Nt Arx- qe do” e (I-B-21) 
Thin 


Under the rather idealized conditions assumed we can obtain a 
formula for stopping power and the statistical variation thereof by 
substituting equation (I-B-14) in equations (I-B-20) and (I-B-21). 
In doing so we may make some assumptions based on the particular 
physical situation we are studying: 

(a) The primary type of particle in the medium responsible for 
stopping effects is the atomic electron. Nuclear collisions do 
contribute slightly to stopping power but such contribution is 
usually small (as will be shown later for our particular case). 

(b) The type of incident particle considered will be heavy, ss our 
work is concerned with protons, 

As a result of these assumptions we may make certain modifica- 
tions in our formulas. If we designate the mass of the electron by 
m , then we may replace Mo by m. The reduced mass m. is also 
approximately equal to m. If we designate the charge of the electron 
by e, then e mu ze and e, is simply e, 


1 


a у 





We see thet the criterion in sub-paragraph (b) on page 10 is 
quite valid. 

If N is defined as the number of m. per unit volume of the 
stopping medium, ve may use it instead of N', provided we also sum 
over all electrons in the atom. Thus N?’ is replaced by N and a 
summation sign over s , where s is the index for esch electron 
of the atom effective in the stopping phenomenon. If there is no 
reason to separate the contributions of various ntomic electrons, or 
if theír properties can be characterized by overall average values, 
then N? can be replaced by NZ, where 7 is the atomic number of 
the element constituting the stopping medium. 


We obtain then for stopping power: 


E = р.н In max 
dx - = T 


nín 
(I-B-22) 
4 
m Y 8 Tain ы 
For stopping power straggling we obtain: 
a 2 

πι: ν΄ 7 μα min? 

TZ DNZT ; (I-B-25) 
max 


provided (as is usual) T »»Τ 


μα» and a is the same for all 


atomic electrons, 
2e Stopping Power - Classical Formulation. 

Inherent in the above analysis is the assumption that the electron 
behaves as a free particle, whereas it is actually bound to the ntom. 
Nevertheless under many circumstances this assumption mav be nuite 


reasonable. For the present discussion let us assume that the perticle 


velocity, v , is greater than the "orbital velocity” of the electrons, 
(The classical approach is essentially the same whether the electron 
is considered as bound by elastic or Coulomb forces to the stom.) If 
u. is the orbital velocity of the sth electron, and I, its ioni- 
zation energy, then 
IF $m u ° : (I-B-24) 
It is possible to divide the stopping power into two ports - 
(a) that portion due to close collisions of the particle with atoms, 
and (b) that portion due to more distant collisions. These must be 
considered sevarately. 
(a) Close collisions are those for which τοι. If Tar? L?’ 
then for most close collisions, the ionization energy is so much 
smaller than T that it may be disregarded. This makes ressonrble the 
assumption that for close collisions the ellos may be considered 
as free. Under such circumstances it is clear that in a direct 
collision, giving an? the eleetron is ejected with a velocity 
practically equal to 2v , so that 
T - Sir (I-B-25) 
Max 
We may note that our original assumption that v > U, leads directly 
to the conclusion that env” >> imu," , Ur T a< > I, - thus 


confirming this assumption. 


Т in for these collisions is of course I, . Then equation 
(I-B-22) leads to 
dE = < 2m т 
[- ST = DN in ~~. ο (I-B-26) 
a 8 8 


(b) For more distant collisions the energy transferred by clsssical 
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collisions would be smaller than the ionization energy. Classical 
models consider the Pes en as an oscillator (real ar virtual) 
which can absorb any amount of energy, however small, when set into 
vibration by a sudden impulse from a changing force field, such as 
would be provided by a passing charged particle. The fact that v 
is greater than u, indicates that the oscillator has moved fron 
its equilibrium position only slightly by the time the particle has 
passed. Thus the "restoring force" has not been mobilized to any 
extent during the collision tims. Once more therefore the electron 
may be considered as undergoing a free collision, under a situstion 
however which permits energy absorption without ejection from the 
atom. 

The value of Tess used ín this region of stop»ing is obviously 


the value of Ta 8 used in the previous case; to wit, І 


4 


The value of T, cannot be zero, since this would give an in- 


in 

finite answer for stopnring power. It has been deemeó necessary to 

presume the existence of an “adiabatic limit" to the distance from 

the atom to the particle path. It 1s presumed that if a particle 

dero at a distance greater than this limit, the electric field is 

applied and released so slowly thet the electron is eased out of end 

back into position without being left in an oscillatory (excited) state. 
The duration of the impulse from the psrticle field is on the 

order of d, /v , where d, is the adiabatic limit corresponding to 


the sth electron. If this duration is less than 1/0 


1 
B» where ω, 


is the oscillatory frequency, then the impulse is faster than the 


ability of the electron to absorb and return it in ediabatic fashion. 


x ^e “ΕΣ 








Thus we make 


ᾱ = Ww, . (I-B-27) 


Then Tain will be the value of energy transfer for collisions at 


this distance. It can be shown that for such distant collisions be- 


tween free particles (N. Bohr, 1948, Eq. 1.1.10): 


κ, 4 
Б 2 z^ e E 1 


T Ва 
ds PUT κ a P ° (I-B-28) 


We see then with the use of equation (I-B-27) that 


° 2 
„ 22% ef Ws г (I-R-29) 


min ἃ 


Equation (I-B-22) leads to: 


4 
Ἢ - DN = ара, (1-B-30) 
dx b 8 2 ze ө“ о 2 


The total stopping power is given by adding the results of 


equations (I-B-26) and (I-B-30): 


-E = 2D0N Ж, merr И (I-B-31) 
8 z e"U& 


This is the classical formule proposed by N. Bohr (1913) . It 
should be noted that the constant k has been introduced arbitrarily 
into the formula. Insofar as the above considerations anply the value 
of k is unity. Bohr'ts more rigorous approach gave to k a value of 
1] 1 . 

3. Elementary Ouartum Considerations in Stopping Power. 

Unfortunately Bohr's formula gave values which were not in very 
good agreement with the experimental data on alpha particle stoning 
which were available at that time. Early attempts to improve the 


formula were unsuccessful. 


---- th кн 





The older quantum theory of atomic structure which Bohr developed 
at about this same time established that the bound energy states of 
the atomic electrons are discrete. This threw doubt upon the classical 
oscillator concept which permits the electron to absorb any amount of 
energy (within the restrictions of the adiabatic limit) and leads to 
formula (I-B-30). Henderson (1922a) suggested that only close , ion- 
izing collisions should be considered, which led to a formula eauiva- 
lent to equation (I-B-26) for the canplete stopping power. This scheme 
likewise failed for alpha werticles, but it is interesting to note that 
it gives answers equal to just about one-half of the experimental 
values (Fowler, 1923) ° 

With our present canprehension of the limitations of classical 
mechanics at atomic distances, it is not difficult to deduce wherein 
formula (I-B-31) fails. We mst accept the quantun requirement that 
energy transfer in a collision rust be at least as large as I, for 
the st electron. (For simplicity, we disregard excitetion possi- 
bilities) Thus classical theory can be expected to hold only if 
T , as given by equation (I-B-29), is greater than I, . This 


min 
eriterion can be transformed by use of the relation 
O T IA (I-B-32) 


and the equation (I-B-24)-into the form 
2 


(3 =- ES . тей 
q κ 


Now if we define x as the velocity of the electron in the ground 
State of the hydrogen atom, we know that 


v = /¥ = 0.219 x 107 em/sec . (I-B-34) 








By substitution in equation (I-B-33) and taking the square root, the 
eriterion becomes: 
ὅτε > v/n ; | (I-B-35) 
ο 8 
Our original assumotion was that v>u, , 50 that all criteria for 


applicability of the classical formula can be expressed as: 





Z Vo Y 
ç > u > 1 . (I-B-36) 


Bohr (1948) showed that this criterion is a little stricter than 
necessary, and that correct results can be expected from the clessical 
formula for values of 2 v/v down to about unity . 

Unfortunately, with respect to protons or alphe varticles, the 
whole stopping power theory herein derived fails when 2 νυν is 
appreciably greater than unity, because the particles tend to lose 
their charge by electron capture under such circumstances (Bohr, 1948; 
also see Section 7, below). Thus the failure of the classical 
equation for such particles is obvious. 

For analysis of stopping power when v >z p (or for protons 
v EN), we must once again divide the effect into two varts. 

(a) Consider the sth electron, for which equations (I-B-24) and 
(I-B-32) are still valid. This electron cannot be vrecisely located, 
but has a odon which is smeared out over a distence having 

"order of magnitude as ° where 

Bp = YU / mu, (1-49 
a, is fairly comparable to atomic dimensions. We first consider the 
contribution to stopping power of those collisions where the minimum 


approach distance of the incident particle is less than a,, which 


g? 
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means that this portion is comparable to part (a) of the classical 
approach. 
For those fast electrons of the stopping material, in case of 


which u, =v, we note that 


particle. ( A equals A/27 , just as M equals h/27 .) 
Since in this discussion the approach distance is less than a, » it 
is always less than A ; and we no longer have a well defined collision 
in the classical sense. Such electrons therefore have little stopping 
effect. 

For electrons in the case of which v>u, (giving Z < a,) we have 
@® situation which still permits a classical apvroach within certain 
limits. 


th electronic cloud 


When the incident particle penetrates the s 
of the atom (that is, roughly, within the distence Ag) the collision 
may be a dírect one, giving Ta = 2 в у. However, many collisions 
will not be quite direct. The least amount of energy able to be 
transferred within this region is indicated by the uncertainty prin- 
cipke. Since we are considering only those collisions where approach 
is closer than a, » our uncertainty in electron position cannot 
exceed this amount. Then equation (I-B-37) indicates that uncertainty 
in momentum cannot be less than m u, ° This in turn means that 


uncertainty in energy is not less than + mu, 2 (=I,)- Thus we can use 


the classical formulas with Tus Pe (Ve are assured that T i, 


A 
— ^ > 
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just es in part (a) of the purely classical approach. ) 


We thus obtain for stopping power in this region a formula which 


is the same as (I-B-26): 


-S2) = on In Zur І-В-38) 
pe ЖИШШ ы: | 


Collisions in this category are called "free collisions" since 
most of them provide energy transfers greatly exceeding the binding 
energy of the electron, permitting it to be considered as free. 

(b) For collision approach distancés equal to P one ean find the 


energy transferred by a formula similar to equation (I-B-28): 


T [L Ë α΄ в“ 1 1 
е5 m v? a Ὁ 


(^ z 3) 
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obtained with the assistance of equations (I-B-24) and (I-B-37). 


(I-B-39) 





This indicates’ that, for heavy particles with velocities well above Vt 
Tas © I, - Thus, if we take a classical point of view, few, if any, 
collisions where the particle passes outside the atom can cause 
ionization. It is this fact which led to the difficulties of 
Henderson (1922) previously mentioned, and such difficulties were 
unresolved until modern wave-mechanics established the possibility 
of energy transfer through the "resonance" process. 

Previous to thís, however, Fermi (1924) had espoused an idea which 
was to be fruitful. He proposed that the perturbing force exerted 
by the particle on the atom could be analyzed as a function of time 


into Fourier components, and the effect of each component on the 


atoms could be compared to the effect of electromagnetic radiation 


-a 





of the same wave-length. ‘This analogy between particle collision 
theory and electromagnetic absorption theory may not be valid for 
close collisions, but for distant collisions where the particle 
field is almost constant over the whole atom the analogy seems 
pronising. Even the idea of the "adiabatic limit" can be carried 
over, since at that remote distance the perturbing field no longer 
conteins oscillatory components of anvreciable magnitude which em 
effect any sort of resonance interaction with the atomic electrons. 
It is well known today that the classical theory of electro- 

magnetic absorption gives results almost the same as those obtained 
by a strict quantum approach and confirmed by experiment (Compton, 
1919; for general discussion see Compton and Allison, 1955). One 
may presume then that for distant collisions we may obtain s correct 
result by continuing the classical approach. This gives us then, 
from equation (I-B-22): 

- &| ° D-N - In δα 


T 


ds 


2 
a 
D:N < fe) , by equations (I-B-28) and 
8 üy 


(I-B-39) ; 
D-N < m ZAE ; (I-B-40) 
8 8 


obtained with the use of equations (I-B-24), (I-B-27), (I-B-32), and 
(I-B-37). 

Collisions leading to thís component of stopving pover are 
Called "resonance collisions", Just as in the case for free collisions, 


those electrons for which u >v ere to be ignored, as their con- 


2 
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tribution according to the formula is quite small. In fact it can 
be seen that the formula even provides a negative (and completely 
unphysical) result when Ig >2 nv? . 


In general, then, the formulation for total stopping power ís 


simply: 
- . rb» —. „an (I-B-41) 
8 8 


being the sum of contributions given in equations (I-B-38) and 
(I-B-40). 

One refinement may be inserted here, for which we shall give 
only a qualitative justification. Hitherto the interaction of the 
atomic electrons has been ignored, but because of the screening 
effects of the electrons on one another they cannot be considered as 
independent oscillators. It is proper to weight the contribution of 
each electron by a factor f » called tie MIEL strength" of 
the s'h electron. The value of each f, is still on the order of 
unity, and it has been shown (Kuhn, 1925; Thomas, 1925; and many 
others) that 

=, & = Z (I-B-42) 

- 
Thus we see that equation (I-B-41) is still approximately correct, 


but more precise answers would be given by: 


. ἀξ E = 2 ву 
| ( I-B-43) 
2 4 
m E S I, 


Except for a slight modification of the term I,, this formule checks 


with that found by Bethe (1930), using s rather strict quantum- 
mechanical approach, similar to that given in the next section, 


е pe. a 
Qt a — 
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We should note here that Bloch (1933) by another rigorous quantum 
mechanical approach obtained a general formula which contained 


Bohr'!s classical eauation (I-B-31) and Bethe's equation (I-B-43) as 
limiting cases, according to whether 2zv ,/V is much greater or 

mich less than unity. (He retained the qualification that w> u, 

for all electrons.) As we have seen, the conditions we are interested 
in are those which put us definitely in the quantum region of appli- 
cability, so that the refinements of Bloch's formula need not be 


considered. 


4. Strict Ouantum-Mechanical Derivation of Stop Power Formula. 





The derivation of stopping power formula using quantum 
mechanics may use any of several different approaches. The one 
outlined here is equivalent to that given by Bethe (1930, 1933). 

The system to be considered consists of an atom and a moving 
charged particle within a large cube of dimension L . The intersction 
between the atom and the moving particle is to be considered as a tíme- 
dependent perturbation. Before the interaction occurs, the particle 
and atom are completely independent, with the former having the 
wave function for a free perticle. The wave-function for the initial 


state of the system can be written therefore as: 


Uy = 5 exp x Po*R ) o » (1-B-44) 


where 
Po = the initial momentum of the "reduced" particle, with re- 
spect to the center-of-mass of the system, = М Y ; 


Н = the reduced mass of the particle-atom system 
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R = the position vector of the particle, referred to the 


atomic nucleus; 


the position vector of the eh eleetron, referred to 


J 
the atomic nucleus; 
Bala = the wave function of the unperturbed atoan in its 
ground state. 

We will concern outselves with the probability of transition 
of this system to a final state in which the atom is left ins 
certain state (excited) characterized by the subscript n (this 
nth state of the atom may be one in which the electrons are still 
bound or may be in that continuum of states where one of the elec- 
trons is emitted fram the atom), and the particle is moving in 
another direction with momentum p. After the collision is 
essentially campleted, the wave-function of the system is egain 
that of atom and free particle without further interaction: 

U, = — exp (ipn) (rj (1-B-45) 

13/2 m - n — 

Let us use the well known formula for first order transition pro- 


babilities (Schiff, 1949, р. 193; Landó, 1951, p. 172) : 


e ET 2 
ET = FA |95.) ; (І-В-46) 
where 
W = the number of particles per unit time scattered into a 


differential element of solid angle (to be precisely 
déscribed below) with final momentum p _, preswming that 
incident particles arrive at such a rate that the re- 


quirement of one incident particle per cube of 





dimension L is meintained; 

РЕ. = the dénsity of final energy states of the system {since 
the atom is left in a definite energy state, the final 
density of states of the system must be obtained by 
consideration of density of states availsble to the 
scattered particle heving momentum p ; the density is 
further curtailed by the fact that we must consider only 
possible directions of the final momentum vector con- 
tained within the element of solid angle specified in 


the definition of WW) ; 


$ 
Hon 


fir U, dt |, Where the integration is carried 
out over all space coordinates of all perticles in the 
system. 

To find % , We note that the number of momentum states avail- 
able to a free particle definitely contained in - cube of dimension 
L, with momentum within the increments dp, ар, япа ер, ís 
(τῷ ιδ] am ар, ара . This ts equally applicable to any element in 
momentum space, so that if we are interested in the number of states 
of the particle with an absolute value of momentum between p= and 
p+dp , the number of momentum stetes available ts (L/h)° p dp 
per unit solid angle. Our solid angle shall be restricted by con- 
sidering only those collisions wherein the scattering anzle in the 
center-of-mass system is betwecn Ө апа Ө t de , so that the solid 
angle is 27T sine: dd. 

The number of energy states in the energy interval dE. is the 


same es the number of momentum states in the momentum interval dp, 


w 





provided the intervals corres»ond. Thus 
Le ЗЕ, = (L/n)Š р? dp (27 sin 8: a6) . (I-B-47) 
The relation between dE. and dp can be obtained from the expression 


for total final energy of the system: 


E Ξ P E tom) 
and 

dE, p 

dp м Y C (1-р-49) 
Therefore, 

3 
Tr Y 
/ E - SCCAR? sin ð de. (I-B-50) 
| hŠ 


Substitution of equation (I-B-50) into equation (I-B-46) gives: 

| M I? # 

i 2 —— — p sine de ES . (I-B-51) 

2 T 
Let us now define d IC) as the differential acattering eross- 
n 

section with angle between Ө апа Ө 4 dO , for one particle inci- 
dent per second upon One square centimeter, and leaving the atan in 
the ntn excited state, This differs from V in that we now have 
l/v particles per cubice centimeter whereas before we had 1/1? 


particles per cubic centimeter. Then 


a $ (9) = (IŠ/v) W > 4 Y 
ο 


2 
- M 15.9 sine ae LA ; (I-B-52) 
оту Po 


This result is the same as that obtained through use of ths Born 
approximation, to which this aprroach is equivalent; and therefore for 


for its validity it depends upon those criteria which make the Rorn 


— 202 {£ ES 
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approximation valid. (Schiff, 1949) 
Let us reduce Hn . The perturbing portion of the Hemiltonian 


is the interaction potential energy between the particle and the atom: 


a 


The first term in the bracket can be neglected, as it gives a result 
for elastic scattering (corresponding to what Bohr calls nuclear 
scattering) which may be ignored for protons or alpha nerticles 
insofer as it affects energy loss, 

It can be shown (Bethe, 1930) that by integration over the 
space coordinates of the particle, Haa can be reduced to the 


following form: 





e? а1 ү 
q? u 2 . € : (I-B-54) 
on 1? |».- p|? n 
where 


28 = JE my [Po - p|: z, ) y. y. E š 


(I-B-55) 


Substitution of the above expression into equation (I-B-52) gives: 


© 2, 
ШИ - STE (É) p emp -a Е. |2 . 
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κ (I-B-56) 
Define a fiew parameter: 
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From the figure we see that 


ge = p^ T pr - 2p Po боз ө. (I-B-57) 
For a given Po апа р, q is a function of 8 only. By dif- 


ferentiation we can derive that: 


sine ag = 29 , (1-B-58) 
P Po 


Then (I-B-56) can be written: 


b Ns 2 
аф (а) = AR. dq БА . (I-B-59) 


This exuression can be written more simply and in a form per- 


mitting easy comnarison with the classical approach given above, if 


we let: 
Po k £ 
2 m om 


Introduction of Q in equation (I-B-59) to replece q as parameter gives: 
3 2 T 2? ef aq ° 
аф (9) = me ae i | | ° (I-B-61) 


We recognize this as being in correspondence with equation 
(I-B-14) obtained by strictly classical means, provided we relate Q 
to T and provided le. is approximately unity for a one-electron 
atom, It is easily shown that Q , fram its definition, is the energy 
the secondary electron would receive from сна collision provided it 
was not bound to the atom, 

To obtein the total probability of a transition which leaves the 


atom in the n™ excited state, we must integrate over all possible 


values of Q : 





X 


2 

pe = з Pr 2 ει . (I-B-62) 

Let us define E as the energy of the atom in its ground state, 
and remember that N equals the number of atoms per cubic centimeter 
of the stopping material. The total energy loss of the particle per 
atom (atomic stopping power) is the average energy gain per atan, 
obtained by multiplying $ by (E, - E,) and then suming (in- 
cluding integrating in the continuum of energy states) over all 
values of n available for excitation. If we multiply by N, we 


then get the total stopping power: 


"max 2 
=. DN (En - E,) 3 i . (I-B-63) 
Smin 
For high values of energy transfer, the ionization energy can 
be ignored; and as we have previously seen, Q Y T. Then for 
iex we may use er Σ ° mw , by equation (I-B-25). (This is 
not accurate for lav values of n since the binding energy of the 
electron has too strong an influence to be ignored in such cases, 
Under normal circumstances these low levels give only a very minor 
contribution to stopping power, so that the aprroximation for Gar 


gives a negligible error overall.) 


For Qi , We cen easily derive an expression from (I-B-60) : 
n 


η 2 
Eun ^ zu πο ο) " 
M , (Ey - Eq) ° ( I-B-64) 


n 4 E 
Further development would be facilitated by reversal of the 
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order of- summation and integration. As the equation stands this is 
not possible since ne is denendent on the particular value of 
the quantum number n . However let us split the integration into 
two perts - with the middle limit, Qo» being independent of n. 
The only other restriction on Qu is that it be less than the 
average ionization potential of the atomic electrons concerned in the 
stopping. I 

Let us determine first of all the sto» pins power formula, 
considering hydrogen as the stonvins medium. If Ry is defined as 
the ionization energy of the hydrogen atom (13.6 ev), we can use 


a =- (Бу)”/ 2m? . (I-B-65) 


Then we obtain: 


= 99, ; e 
“ах a oe = (E - Eo) [€ | 
lo 
« 2 
É a gc) 2 [ει] . (1-5-66) 
Amin 


In the evaluation of the first integral, we can use a theorem 
proved by Bethe (1930), which indicates that = (E, - E) ll a 

Evaluation of the second integral can be simplified since 0. and 
чыл are both much less than Ry. This indicates that © is very 
small; condequently p,- p is small and approximately perpendicular 
to the path of the incident particle. The situation becomes closely 
equivelent to the case of energy absorption from electromagnetic waves. 
In auch cases the square of the matrix element E, can be shorm to 


equal the optical oscillstor strength for the transition from ground 


-- 20 = 





to n"! excited state of the atom. 


Bethe (1930; 1933) uses this 


simplifying fact and carries out a complete numerical solution for 


the stopping power of hydrogen. 


In doing so he uses the well known 


wave functions for hydrogen - expressed in polar spherical coordinates 


for the discrete levels, and in parabolic coordimtes for Bhose levels 


He obtains: 


Q 2 


in the continuum above. 
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= .- 0:098 . (I-B-66) 
- 0,098 
2 
2nv . (I-B-67) 
1.103 Ry 


‘This formula is fairly accurate and is eonfirmed avproximately 


by experimental work (Mano, 1934). 


It should be noted that in its 


derivation it has been necessary to do the same thing ss vas done in 


sub-sections 2 and 3 above - that is, to divide the stop»ing effect 


into two parts, that involving close collísions where the binding 


energy of the electron is negligible, and that involving distant 


collisions where anslogy to electromagnetic dispersion and absorption 


theory is utilized. 


For heavier elements the situation is a little more complex. 


Bethe (1930) shows thet for each electron of the stopping m*terial, 


the above formula is valid, provided that the factor 1.103 Ry is 


replaced by a term I: » which is somewhat larger than the ionization 


energy of the electron. 


— 724 = 


Further, in suming the effect of each 





electron, it must be weighted by its oscillator strength. 
The determination of I! is rather complex, as its precise value 
depends upon a knowledge of the wave-functions for the etom concerned. 


The formula for complex atoms may be stated as follows: 


dE 477 22 el y 
Eo. fT = t, n2BY | (1-5-68) 
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Тв 
4, 
- 4T z^ e* nu Z7 ln £r Y : (I-B-69) 
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where I ç » called the "average excitation energy," is theoretically 
defined by: 


zul, = E tdg» . (I-B-70) 


Formula (I-B-68) is of course identical with (I-B-43) above, except 


that I* is probably a little larger than I the true ionization 


8 g» 


energy. In practice it is necessary to select empirically the value 
of I y or the values of one or more IS to make the formula agree with 
experimental data. 
9. Modifications to Theory for Slow Particles Through Heavy Media. 
Most of the above theory was developed on the basis that the 
eharged particle has a velocity faster than the "orbital velocity" 
of each of the atomic electrons in the stopping medium.  Th»t the 
criterion fails for protons having energies íin the region we are 
studying is easy to show. If we consider the velocity of an atomic 
electron given by equation (I-B-24), then the energy of a proton 
having the same velocity as an atomic electron of 1onization energy 
I, is given by (11/m) I, . he eriterion for validity of the unmodi- 


fied equation for stopping power becomes: 
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E >> (M/m) 1, for every electron of the atom, where E is the 
actual proton energy. 
From tables of idnization energies (Siegbahn, 1931), it is easily 


computed that the criterion fails: 


For „N , K electrons, when Е Ж 0.72 Mev; 
For go ‚ К electrons, when E C 1.0 Mev; 
For 1018 K electrons, when E = 1.6 Mev; 
For A „K electrons, when E Ξ 5.8 Mev; 


18 
For 194 » K and L electrons, when E Z 0.35 Mev; 


For >Ni, K and L electrons, when E = 1.5 Mev; 


F Kr, K, L, and M electrons, when E Z 0.48 Lev; 


Or ag 

For g,Xe, K, L, and M electrons, when E % 1.5 Key, 

What does one do 1f the criterion fails? One approximation 
(Livingston and Bethe, 1937) involves the assumption that all elec- 
trons for which E > (M/m) I, are fully effective in stopping, and 
the other electrons do not contribute at all. This means that 


formula (I-B-69) is usable provided 2' is used in place of Z, 


where 


reb m To, (1-B-71) 
n 


the sumation being taken only over those electrons slower than the 
incident particle. -— is also different in such cases, since only 
those outer electrons conforming to the criterion are used in the 
averaging. | 

This method gives a curve of proper shape but requires experimen- 
tal values at various energies to provide precise accuracy, Further- 


more, unless Lay is readjusted accordingly, discontinuities are likely 





to appear in the curve at the points where T = (N/m)I, , since 
Z* changes at these points. 

Bohr (1948), who was less interested in precise resubts than 
in general trends, evaded this difficulty. He recognized in his 
analysis which led to the eauation called herein equation (I-B-41) 
that those electrons for which u, >Y should be ignored. To 
determine those electrons which are effective in the stopping process, 
he used a formula based on the Thomas-Fermi statistical model of the 
atom (Thomas, 1927; Fermi, 1928): 

n(u,) = the number of electrons in the atom with velocity 


smaller than Ug, ; 


πὲ 


215 ῃ/π. . (1-B-72) 
Equation (I-B-41) is modified by integrating over dn(u, ) instead 


of summing over s : 
2v 
. 2E 2DN A àn(u) . (1-B-75) 
sau, (= lg) 





0 
The upper limit is that value for which the logarithm goes to zero. 


Simplificetion and integration then gives: 


8px zY5 (v/v,) 


dx 
_ 16722 tN ae (I-B-74) 
= m v, σ 


Bohr states that the formula he uses for n(u,), is not very accurate 
if u, represents the velocity of the very outer or very inner elec- 
trons. Thus formula (I-B-74) may be expected to be most valid for 

those particle velocities which are comparable to electron velocities 


in the intermediate shells. 


Neufeld (1950) attampted to improve Bohr’'s annroach во ас to 
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provide answers more accurate tnan that given above. He went back 
to the basic formuletion wherein stopping effects are separated into 
two parts - the contribution from "free collisions" ond that from 
"resonance collisions". For the latter he used a more elaborste 
formula which had first been proposed by Fermi (1940). This formula 
is noteworthy since it provides positive contributions from every 
electron, regardless of velocity. Furthermore it is convergent if 
all distant collisions are included, even out to a collision distance 
of infinity. This obviates the need for the "adiabatic limit" 
concept. 

Neufeld'!s approach is tedious and requires a separate set of 
calculations for each type of stopping element. His calculated 
values for a few metals seem to be about 10% less then the best 
experimental results available at the moment. The advantage of his 
approach is that it recuires no previous experimental work to 
determine the parameters. 

We have reserved till last the discussion of a more exact 
treatment of this problem, which was initiated by Dmm end 
Bethe (1937). As a prelude to this contribution, it should be noted 
that shortly after Bethe (1930) gave his stopping pover formula, Mott 
(1931) and Henneberg (1933) showed that the Born method of treating 
collision problems could be applied with validity to collisions of 
heavy particles with heavy atoms having electronic velocities foster 
than the mrticle velocity. This is true for atoms of high atomic 
number, since the perturbing potential of the incident nssticle is 


small compared to the binding potential for the fast, inner electrons. 
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(See also Nott and Massey, 1949.) This fact gives one faith in the 
use of previously derived formulas even when the original criterion 
that v> u, is not completely adhered to. 

Suppose one wishes to obtain results which are valid for par- 
ticle velocities dovm to the order of magnitude of the velocity of 
the K-shell electrons, or somewhat lower, An exact solution of the 
contribution of the K electrons is obviously much better than an 
assumption requiring either full contribution or zero contribution 
of the K electrons to the stopping effect. e may note that fairly 
Simple wave-functions are available for use here, since even ina 
complex atom the innermost electrons obey wave-functions which are 
almost hydrogen-like, The formula used by Livingston and Bethe is 
thus simply a modification of the equation used by Bethe to derive 


the stopping power of hydrogen (equation I-B-63): 


4 
|- ΕΙ pS Rex S , (I-B-75) 
where | кы Е 
B. = Er an: e.g. (1-B-76) 
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In this formula one integrates instead of summing over the 
various energy changes, since most of the discrete levels well below 
the continuum can be considered as filled. The scele of the values 
of energy and Q are changed to put all calculations on a uniform 
basis for different values of 2: 


E -E 
Et - ET Wt RS ; (I-B-77) 








where Z off = effective nuclear charge as felt by the K electrons, 
= Z - (inner screening constant) (Slater, 1930) 
E — 
Ὁ- = 2 ; (I-B-78) 
Zeer Ὁ 5 
(Po - P) 
Q cz — : BD ; (I-B-79) 
Lett eRY on 
fe: E 
ζ = 2 Т = : , (I-B-80) 
Zeer CRY eff “Ry 
mins Bean. (1-B-81) 


which is actually the same ss eauation (I-B-64). 

@ is the lower limit of E*' , inasmuch ss we consider transitions 
to the discrete levels as impossible by reason of their being filled. 
Infinity is permitted as upper limit for Q, since the computations 
of the matrix element make the probabilities exceedingly small for 
any contributions by values of Q above env”, 

Livingston and Bethe have evaluated this formula for B. by 
numerical integration. For © equal to 0.7 , which is epnroximately 
correct for values of Z from 6 to 13, they have plotted By against 
n as argument. B, can also be expressed as: 

By, EE ^) = ап 3.63N - οκίή le (I-B-82)7 
where C. is a corrective term to the basic formula (I-B-68 or 
I-B-69) whieh approsches zero as particle energy (or ? ) increnses. 


Curves for C are also given by Livingston and Bethe. 


н It is easy to show that the factor 1.81 is approximately the sum 
of the oscillator strengths of the two K electrons, while 
3.630 Z u-—-— Be . : 
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They have also given a curve for atomic stopping power of air for 
protons, data from which is tabulated in Table I, over the energy 
range in which we are interested. 

This formulation of K electron stopping power has been ex- 
tended by Brown (1950) and l'alske (1952) to values of $ higher 
than 0,7 , which in effect extends the theory to heavier atoms, 

They have made the computations and supplied graphical results for 
Bk( € ,4) and e Co) ° 

It is pertinent to ask about the corrections required for 
L, M, and higher electrons as the velocity of the incident particle 
approaches their orbital velocities. This is porticulerly impor- 
tant, as we have seen, for the heavier elements. The exact analysis 
for stopping power of electrons in these shells is well nigh im- 
possible to carry out because of the complexity and lack of knowledge 
of ihe vave-functions ot these shells, 

Hirschfelder and Magee (1948) generalized the foregoing method 
to apply to shells higher than the K shall, ‘They utilized the 
original curves provided by Livingston and Bethe since those of 
Brown or VWalske were not availeble at that time, In effect they 


notec that tne stopping number B , where 


a = 2 in 2a Cy - (higher shell corrections), 
av (I-B-83) 
can be expressed as: 
"es Murat, ϱ (I-B-84) 
1 "1 1 
where 
= number of electrons in the th shell; 
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τι = effective oscillator strength of an electron of the 
ith shell; 
Pa = the essential variable for the 1 shell, depending 


upon N : 


For the K shell it is seen from enuation (I-B-82) that 


BAR ) Cr 
Brus τι ο της ° ΡΝ 
and, as before, 
E el : (I-B-86) 
l  (M/m) Iw 
K 
with 2 
Ір = Zee Y ° (I-B-87) 


: The generalization was made in the assumption that, in forn, 
δ, is the same as be for'ell values of 1. It is only necessary 
to use a different value of rt (the ionization energy without 
external screening for the shell considered) in computing the 


argument ? with which to use the curves for B, . For all but the 


outer shell, the value of I'* was taken to be: 


Y Y = Me Ἢ ( I-B-88) 


with n4 being the effective principal quantum number of the ¿th 
shell. 

It is hardly necessary to point out that such a procedure is 
approximate - primarily because of the marked deviation of the wave- 
functions of the middle shells from a hydrogen-like character. The 
originators of this approach hoped to predict the trends rather 
closely however; and the semi-empirical character of Bethe's original 


formulation was retained by adjusting the value of I'' for the outer 
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Table Τ. 


Semi-empirical Calculations of Atomic Stopping Power 


Air - Livingston and Bethe (1937) 
Argon- Hirschfelder and Magee (1948) 
Xenon- Hirschfelder and Magee (1948) 







O (10715 ev-cm2 ) 


Air Argon Xenon 









Dos 12.95 20,834 35.5828 
0,4 10.63 17.601 30.6465 
aÑo 8.86 15.432 27.2376 












0.6 7.70 13.853 24,7589 
0.7 6.90 12,604 22,7846 
0.8 6.30 11,610 21.1980 
0.9 5.81 10.792 19.8589 
1.0 9.45 10.108 18.7458 







27922 17.7874 
4.80 9,011 16. 9514 
4.55 8.563 16.1795 
4.29 8.160 15,5181 
7.807 14,9227 
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shell of electrons as necessary to conform to experimentel values of 
range in the media (gaseous) selected for study. 

The studies of Hirschfelder and Magee were applied to two gases 
of particular interest here - argon and xenon. The computed results 
for atomic stopping power of these two gases for protons are tabu- 
lated in Table I , over the range of enim en in which we are interested, 
6. Effect of the Physical anà Chemical State of Vedia on Stopping Power. 

In all of the theory presented previously the atom which absorbs 
the energy has been considered to be in a free state, such as a 
monatomic gas. Any sort of chemical combinetion or condensation of 
the medium into a liquid or solid state might be expected to shift the 
energy levels, at least of the electrons in the outer shells, and thus 
affect the stopping power by changing the value of Iv , the average 
excitation energy. 

As far as chemical combination in gases is concerned, the effect, 
if present, must be auite small. If the atoms act independently, the 
stopping power per molecule should be a simple sum of the atomic 
stopping powers of the individual atoms. This rule seens valid in 
gases though not always in linuids and solids. (See the excellent 
reviews of this aspect of stopping power given by Geiger (1927) and 
Taylor (1952).) We are concerned with the stopping power of com- 
pounds in two respects. First, some of the impurities found in the 
gases are compounded of more than one element; second, one of the 


primary gases studied in this research (nitrogen) occurs іп а diatomic 


The relation between range and stopving power is easily seen to be 


Eo an 
R(E) = - (I-B-89) 
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state. 

Theory hints that more appreciable effects might be forth- 
coning if the stopping medium is in a condensed state, because of 
the influence of the many adjacent atoms on the outer electron energy 
levels. Some early experiments seemed to indicate a consistent 
difference between the stopping of alphs particles by elemental 
solids and gases (Brizgs, 1927; Geiger, 1927). These seeming 
differences were later shown to be due to experimental factors 
(Livingston and Bethe, 1937). Experimental data - not entirely in 
agreement among themselves - still accumulate however to indicate 
that special effects may be present in certain coses (Taylor, 1952). 

Both theory and experiment khsve shown irked effects in 
condensed substances at relativistic velocities of incident n-rticles, 
but we sre not interested in that energy region. 

Theoreticians have given special attention to stopping by metals. 
The outer electrons in metals are practically free and their behavior 
is radically different from those electrons bound to individuel atoms. 
As early as 1933, von Weisacker (1933) derived a formula which strongly 
related stopping power to conductivity. His theory was shovn by 
subsequent experiment (Gerritsen, 1946) to nake incorrect predictions, 
and Kramers (1947) proved theoretically that special conduction 
effects were practically cancelled by polarization effects, A. Bohr 
(1948) and Pines (1952) published further estínates of small reduc- 
tions іп stopping power by free valence electrons for incident particles 
of moderate velocity. These corrective factors are not large, however, 


and for heavier elements wherein the proportion of conduction elec- 





trons is small the corrections predicted are auite negligible. 

This is confirmed experimentally by Heller and Tendam (1951), who 
found that for Z213 no consistent differences between stopping 
power of metals and semi-conductors are found except the expected 
variations with atomic number. (Note should be made that such 
results were obtained only for perticle velocities in excess of the 
velocity of 1 Mev protons.) The only experimental confirmation of 
special effects due to conduction properties of metals is the work 
of Madsen and Venkateswarlu (1948a) on sto»nving by beryllium for 
protons in the energy region from 500 to 1500 kev. The effect is 
expressed as an abnormally lerge value of I v for beryllium to be 
used in the stopping power equation. 

A. Bohr further indicated possible polarization effects mey 
affect stopping power for condensed media even when non-conducting. 
Halpern and Hall (1948) also suggested the possibility that the more 
loosely bound electrons in condensed media may have different oscil- 
lator frecuencies associated with them than hhey would have in the 
gaseous state. The special effects predicted are still small except 
in those cases where the outer electrons constitute a large pro- 


portion of all electrons effective in stopping. 
7. Capture and Loss of Electrons. 

In 1922, Henderson (1922b) observed that a beam of alpha porticles 
contains a portion of singly charged particles, which portion increases 
with decreasing beam velocity. Since thet time experiment and theory 
have confirmed that for all charged particles there is a pronounced 


tendency at low velocities for the particle to change its charge by 
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capture and loss of electrons. It is desirable to consider the 
phenomenon briefly in respect to stopping power of protons, since 
the stopping power theory presented above depends on the retention 
of positive charge by the proton. 

Theory (Bohr, 1948) has given an approximate answer to the 
cross-section for capture of an electron by a proton, giving neutral 
hydrogen; and also a formula is provided for the cross-section for 
loss of this electron by the neutral hydrogen in flight. A proton 
undergoes in its passage through matter a succession of neutralizing 
and re-ionizing events, so that the proportion of its flight which 
occurs in the charged state may be obtained by taking the ratio of 
eross-sections. 


Let y = the proportion of flight as neutral hydrogen; 
Te 


1-Β-90 
σε - ση , ( ) 
where 
σ = һе electron capture cross-section ; 
ο. = the electron loss cross-section . 


If 01 >> O% » then 
y 2 9/9 . (I-P-91) 


According to Bohr (1948), 


Y = AT a 7 2° 21/3 (v/v) ® : (I-P-92) 
E. 2 ¿-1 „2/3 
ES = Te, $c Z (v/v) Е (I-R-93) 
Since z equals unity, 
у = 4271/5 (v/v)? . (I-B-94) 


Now in our experimental work the lowest value of v is 0,81 x 10? 


cm/sec (corresponding to an energy of 340 kev), the lowest value of Z 
— < - 





is 7, and the value of " is 0.219 x 109 em/sec , so that Lt equals 
0,0031. Thus y is so low and decreases so repidly as particle 
velocity increases that we may ignore entirely the possibility of 
electron capture. 

A confirmatory note should be added. T, Hall (see Warshaw, 1949) 
has found that experimental results are in substantial agreement with 
Bohr*ts theory used above, 

B. Straggling in Stopping Power. 
Bohr's equation for straggling in stopping power is given by 


substituting equation (I-B-25) into equation (I-B-23): 


(N = 417522 uz. (І-В-95) 

This eouation is valid for fast incident particles. For slower 
particles, such that some atomie electrons sre no longer slow compared 
to the particle velocity, Bohr returns to the Thomas-Fermi model and 
uses only those electrons slower than 2v . That is, Z is replaced 
by Z*, where 

2 = 21/5 (2v/v,) > (I-B-86) 
as determined from equation (I-B-72). 

Livingston and Bethe (1937) outline an apvroach to this prob- 
lem which parallels the analysis of Bobr, previously given, in much the 
same way that Bethe's derivation of stopping pover formula is analo- 


gous to Bohrts classical apnroach. They obtain: 


q е 
ZÉ) = 4% 2% of a+ SK In ^a ITA 
x n f 
n mv? I 


(I-B-97) 


where Z* is defined ass in equation (I-B-71); Kn is a numerical 
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factor usually equal to about 4/3 ; 2, is the number of electrons 
in the neh shell; and In is the average excitation energy of the 
nth shell, as usual. 

This formula is somewhat difficult to use, as the values of In 
can only be estimated. Furthermore the value of ¿E given ís only 
approximate and its true value is difficult to determine. 

The extra factor in the bracket is the essential correction to 
Bohr's formula, and actually seems to be not so much a true quantum 
correction, as a term to to take into account the motion of the 
orbital electrons. (The factor i in the numerator is actually 
only an aporoximation for the kinetic energy of the electron.) 
Bohr*s theory, it must be remembered, assumes that the electrons 
are motionless before collision. Averaged over all particle- 
electron collisions, this is approximately true, so that the stopping 
power relation is not affected by such an assumption; but energy 
tSansfer variations in individual collision cases are enlarged be- 


cause of this electron motion, so that the stragrling is increased. 


C. Review of Existing Experimental Data. 
1. Stopping Power. 
To review all experimental date which has ever been published 
on stopping power would be very tedious and quite unnecessary. At 
the opposite extreme, to list available experimental results only if 
they cover precisely the same set of conditions and materials which 
we ere studying would be rather simple - there is very little to date. 
It is profitable to take a middle course - to list not only that 


data which may duplicate our own, but also to list data in either one 
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of the following categories: 
(a) stopping power for protons of the materials we are study- 
ing in energy regions to either side of the regions we are 
concerned with (400 ~ 1100 kev); 
(b) stopping power for protons of other elements, in the energy 
region which coincides with ours. 

In this way we may synthesize our results with others into the 


overall stopping power scheme. 


Table II. 


Experimentally Determined Values of 


Average Excitation Energy (I) 





Element! 2 | Zt (if Particle and 
Е а velocity ог Reference 
energy 
sN, 7 81 * lano, 1934 
Ne 10 132 ο и " 
Α 18 195 5 " п 
Ni 28 325 * N " 
Cu 29 320 * " " 
279 300 Mev prot. Bakker & Segré, 1951 
307 240 " " Mether & Segró, 1951 
375 7ο m" z Bloembergen & van 
365 100 " E Heerden, 1951 
27.4 293 Py" " Liv. & nethe, 1937 
Kr 36 390 * Hano, 1934 
Xe 54 530 * " " 


Alpha particles from natursl radioactive elements (6 - 11 Mev), 
which corresponds to protons having velocities equivalent to 
1.5 - 2.8 Mev of energy. 


— < 


Let us first list all data available for energies higher than 
the range studied herein. This is best summarized by indicating the 
value of I ç found necessary to make Bethe's formula (equation I-B-69) 
cheek with experiment (see Table II, above). 

On the low energy side of the energy spectrum, some data is 
available relating to the materials in which we are interested. A 
certain amount of work hes been done by the "Kevatron" group at 
Chicago. Wilcox (1948) published results at low proton energies for 
aluminum and gold, but his results have been superseded by better data. 
Warshaw (1949) has published data for protons of energies from 40 to 
360 kev for aluminum, copper, silver, end gold. He extrapolated his 
curves up to several Mev, where the theoretical formule of Bethe would 
be expected to be valid. 

Work for gases has been reported recently both from the Univer- 
sity of Chicago and the California Institute of Technology. From the 
latter institution Dunbar, et al. (1952) have reported on stopping 
power of some gases for vrotons of energies from 25 to 600 kev. 
Additional data has been recently provided by this group of workers! 
which means that most of the simple gases, ineluding all those we 
are reporting on herein, have been covered for protons up to 600 kev 
of energy. Their curves sre reproduced in Figure 13, to show the 
match with our data. Weyl (1953) at Chicago hes provided data for 
stopping power of various light gases, including argon, for protons of 


energies from 50 to 500 kev. His data is consistent with that from 


l up-to-date information from the Cal.Tech. laboratory has been kindly 


transmitted by Prof. WJ. Whaling (personal communication). This is 
expected to be published later under authorship of Dunbar et al. 
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the California Institute of Technology and ds therefore not plotted. 1 

Stopning data providing a complete coverage of the particular 
region of interest to us has been provided by only a few workers. 
Madsen and Venkateswarlu (1948a), being the first to employ the 
experimental technique we employ with the Van de Grasff generator, 
published results for beryllium stopping of protons with energies 
from 500 to 1500 kev. Using Bethe's formula, and including the 
correction term E they found the data to be consistent with a 
value of I ç of 64£5 ev, 

Using the same technique Huus and Madsen (1949) obtained a 
little data in thie region for gold. Madsen has also recently sub- 
mitted some data for copper.” His results for the latter element are 
not consistent with other experimental work, which is possibly due to 
difficulties in obtaining foils of consistent thickness. 

The greatest amount of work in this region is provided by 
Kahn (1953), who has obtained curves for protons of energies from 
400 kev to 2.0 Mev for beryllium, aluminum, copper, and gold. 

Figure 12 contains his results, which agrees with the earlier data on 
beryllium and gold mentioned above. 
2. Shraccling. 

Very little previous experimental work has been published on 

straggling in stopping power in the energy region of interest to us. 


The reason is easily seen. Even the techniques for stopying power 
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f Detailed information on the work ət Chicago, es well as other 


information available to him, has been kindly provided by Pfof. 
S. K. Allison (personal communication). 


? ue acknowledge appreciatively a personal communication from Dr. 
Madsen, giving detailed information on all his work to date. 
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determination in this energy region are somewhat crude, and exneri- 
mental variations sometimes wide. It can be appreciated that the 
theoretically valid deviations in stopping power can be easily masked 
by His Sumo factors which cause inaccuracies in stopping vower ex- 
periments. 

She only experimental work published on straggling is by Madsen 
and Venkateswarlu (1948b), who were able to analyze the results of 
their stopping power data (1948a) to give straggling in adition to 
mean stopping power. Results are for beryllium, with proton energies 
from 400 to 1200 kev. The results vary widely, but seen on the average 
to be in reasonable agreement with the theory of Bohr (1948), given 
herein as equation (I-B-95). The value of Dan? was found to 
average 8,99 De? , with individuel resdings ranging from 6.9 


to 10.4 . In these units, the theoretical prediction is 8.3. 
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TI. Basic Experimental Technique. 


A. Basic Principle. 
If a beam of monoenergetic protons impinges upon a target of 


a certain material, a nuclear reaction may be induced which pro- 
duces quante or perticles as a result. The investigation of proton- 
gamma reactions requires the determination of the number of garma 
photons produced (and detected) per designated number of protons 
incident on the target material. For a given target, such a figure, 
which we may call the yield, is a function of the energy of the in- 
eident proton. By varying the proton energy in small increments and 
messuring the yield at each step, one can accumulate data on yield 
versus proton energy. For meny light isotopes as targets, a vlot of 
these results in the Van de Graaff energy region (a few hundred kev 
up to several Mev) gives a type of curve which is low end slowly 
rising, on which is superposed a number of "resonance Peaks" whose 
height, width, and position in the energy spectrum sre charscteristie 
of the isotope(s) involved. 

How if a given thickness of matter is interposed in the vath of 
the beam after it leaves the accelerator but before it hits the target, 
the protons are slowed by amounts, the average of which is determined 
by the average energy loss of the protons in the interposed material. 
Then to duplicate a "peak" in a yield-versus-proton energy experiment, 
the protons obviously must be given an initial energy greater than that 
characteristic of the target material by just the amount of energy loss 


in the interposed material. 
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If, in such ar experiment, one measures proton energy before 
the beam strikes the interposed stopping material, a yield-enerry 
curve may be constructed, called hereafter a shifted or displaced 
epectrum. The peak-to-peak correspondence with the undeviated spectrum 
is easily established by inspection. 

The undeviated peak energy values may be determined bv exveri- 
ments without stopping media interposed or, if already well estab- 
lished, they may be determined from the literature. One then can 
subtrect the energy values for the corresponcing peaks in the 
deviated and undeviated spectrum to obtain energy loss AF in the 
material. Dividing by the path length of the protons in the medium 
gives AE/Ax . The theorem of mean velue indicates that this v-lue, 
which depends on the energy of the beam as it enters and as it leaves 
the medium, is the same as the velue of dE/dx for some value of E 
between the incident and exit beam energy. This mean value of E is 
not difficult to determine when AE is small comvared to E. 

It is found that the peaks in the displaced spectrum are brooder 
thon those in the basic, unshifted spectrum. This is cnused by 
individual proton variation, or straggling, in stonning power. hat 
is, when deviations in energy loss are possible for individual protons 
passing through matter, then beams of averege energy anrreciably 
different from the proper displeced peak value will still contain 
protons of the proper energy to excite the target atoms in the 
resonence region. 

The technique briefly outlinedabove was first used by Mndsen and 


Venkateswarlu (1948a,b). A more detailed analvsis is now presented, 





showins the corrections necessary to get = true answer from experi- 
mental results, 


Β. Detailed Analysis. 
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Figure 2 


Figure 3 is a shhematic diagram o? the annaratus for measure- 
ment of stopping power, with the voltages specified for a proton 
fron the moment it leaves the accelerator section until it enters 
the target nucleus to release n phóton to be counted. In this 


figure, 


E energy of a photon-producing proton as it leaves the 


O 


accelerator; 

E, = energy of the photon-producing proton efter nessin- 
through the analyzing field; 

E = energy of the photon-producing proton after pascing the 
defining slit and as it enters the foil; 

E z energy of the photon-producing proton ss it leaves the 
foil and enters the gas spoce; 


energy of the photon-produeing »roton as it leaves the 


{51 
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gas space and enters the target moterial; 


= 


- 
£ 
~ чь 
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Ἐς = energy of the photon-producing proton as it enters the 
nucleus to produce the reaction. 
le note that if the gas chamber is evacuated, Ey = Pa ° 
In order to find the energy E¿ of a marma-producing proton as 
it comes fron the accelerator, it is possible to start with Es and 
successively add to it tho energy losses suffered ín all verts of the 


system between E. and E, « Thus 


5 | 
E - EgdAx|«^ E ΑΣ ΙΔ E 5251. (II-B-1) 


where 
Alp = the energy change of the proton in passing through a 


certain amount of the target material; 
AE = the energy change of the proton in passing through 
the gas; 


the energy change of the proton in passing through 


- 
the foil; 

AE. ms the energy change of the ppoton in passing through the 
defining slit, or as a result of the regulatory action 
of the defining slit; 

AE, = the energy change of the proton in passing through the 
magnet, or as a result of any miscellaneous effect in 
the systen. 

What we are particularly concerned with, however, are the 
relationships among the distributions in energy of the probability of 
detecting a garma photon for each proton accepted from the accelerator. 
These distribution functions depend not only upon proton energy, but 


also upon the position in the apparatus where the energy is specified, 
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so that we conclude that there is a probability distribution for each 
of the variables in equation (II-B-l). 

The statistical variables of interest are the mean and the 
variance (square of the standard deviation) of the distributions. For 
the mean value: 


Em Es + [Ap + [au | + las; + [ae |+[45, | ; (II-B-2) 


and for the variances, the following formula is valid (Hoel, 1947) : 


2 2 2 2 P 2 2 
N, i Le +0, UNE K Ἢ +A, , (II-B-3) 
where the subscripts refer to the same portion of the apparatus as 

in equation (II-B-1). The precise values of some of these statis- 
tical variables will be determined later, but others can be simplified 
in such an obvious way that we may conviently do so at this point. 

We assume that the resonance peaks are so widely separated that 
they do not interfere with one another, and thus thé statistical 
analysis may be carried out on each peak alone. Furthermore, the 
frequency distribution is considered sufficiently symmetrical that the 
maximum value can te taken as the mean value without appreciable error. 
То simplify the notation let us assume, unless otherwise stated, that 
when we refer to energy henceforth, we mean the mean value of the 
energy of the whole proton bean, without the necessity of placement of 
a bar over the symbol. 

Eo becomes the average beam energy from the accelerator giving 
the resonance maximum. Es is the value of proton energy as it enters 


the target nucleus to give the resonance maximum, and we will therefore 


re-designate it as Bnat’ AEn ís the average value of energy loss of 
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protons in the target material, and is therefore #7 , where T is 
the target thickness expressed es the average energy loss for 
protons passing through the target. If we assume that the effect 
on individuel protons b» the slit and the magnet sre as likely to be 
in a vositive as a negative direction (see discussion below) then 
AE, and ДЕ, аге zero, | 

1 


resonance maximum point, and is determined fron the experimental data; 


ο ig the standard deviation of the yield curve about the 
whereas ο is the standard deviation of the natural, or theoretical, 
curve and will be called e henceforth. N, and TE. mensure 
the straggling effect in energy loss ín the gas and foil respectively. 
The deviations caused by the slit and magnet may be exvlained somerhat: 
(a) The slit is of finite width, and the voltage regulation system 

of the Van de Graaff generator (Smith, 1952) is such as to keep the 
proton beam between the sides of the slit. Then the maximum swing 

in beam voltage is on the order of the amount required to swing the 
beam from one side of the slit to the other. If the value of average 
beam voltace is measured with the bean centered in the slit, the 
deviations ere equally likely to be plus or minus. 

(b) We may expect other miscellaneous variations of s presumed random 
nature: inhomogeneities in the magnetic field, variations in the 
vertical end horizontal position of the princival axis of the focused 
beam, slight variations in the electromagnetic field current, and 

ШШ | ›лпесшв variations of a human or instrumental nature. These are 


eonveniently grouped together with the subscript M ° 


Using the modifications given above we can now re-write equations 


= 


€, 
_ u 


<= 





(II-B-2) and (II-B-3). Let us first do so for an experimental 
curve without any stop»ing media at all: 


BL re Et Ar ; (II-B-4) 


nd 2 2 2 2 2 
= E) + + ὃ 
iN 0 nat T 0 5 N wo (II-B-5) 


Then we can write the equations valid for experimental curves where 


a 


stopping media are interposed in the beam. In doing so those values 
which are the same for the corresponding peak without media inter- 


posed are left unvrimed, values which are shifted are primed: 


aa = fo, + + [^x.| *] ^5. | ; (II-B-6) 
and 
‚2 2 2 2 ‚2 ‚2 
= + + + + БЕ 
I, A í 4 (1 Ir Ts (II-B-7) 


In the latter formula, Te signifies the sum of the variations in 
energy loss througk both foil and gas. In practice it will denote 
only the variation for the foil, inssmuch as we have obtained no 
results in cases where both gas and foil are used, 

Subtracting equations (II-B-4) from (II-B-6), and (II-B-5) from 
(II-B-7) gives the basic equations which permit determination of 
stopping power and straggling from experimental results for the 
resonance peaks: 


Without gas, 


3 - . nc 
A, = ES ES ; (II-B-8) 
with gas and foil, 
- 2? = = 9 - Ң = 
AE, = ES ES AE, ; (II-B-9) 
and without gas, 
2 ‚2 2 ‚2 2 ee 2 
N = (fl, - £1) - (05 - Ως) -αἰμ = An 
(II-B-10) 


That the primed and unprimed values of ΠῚ g and τα do not always 
cancel may be understood from the discussion en sections III.F.3 and 4, 


-bf == 


111. Experimental Details and Preliminary Calculations. 


The primeiple of operation of the Van de Graaff type of 
generator-accelerator is too well known to require any explanation 
here. A description of the particular one at the Ohio State 
University and its associated auxiliary equipment may be found 
elsewbere (Grove, 1947; Cooper, 1949; Grove, 1950; Taylor, 1952; 
Smith, 1952). Certain improvements and modifications wére 
required in order to fit the machine to the purpose of the present 
researeh. They will be described in the sections immediately following. 
А. Target-Faraday Cage Assembly. 

In order that the protons may be collected and counted, the 
target on which the protons impinge is mounted on an assembly called 
a "Faraday cage." This is insulated electrically from the rest of 
the accelerator, so that the collected positive charge can drain to 
ground through an "integrator" which measures this charge. 

It was necessary to provide a mechanism for introducing a foil 
into the proton besm; also a gas~-holding chamber was required with a 
thin foil-covered window so that protons could be introduced to 
traverse the gas. For this purpose a tyne of valve box” ms 
modified to perform all the functions of target holder, Faraday cage, 
gas holder, and foil retainer. This assembly is depicted in Figure 
4, and is almost self~explanatory. 

The system may operate with the valve down for basic calibration 


runs. If stopping power of a foil is to be investigated, the foil 





E This box was originally made after a type used on the University of 
Wisconsin Van de Graaff generator ghñe pattern for which was 
generously provided by Prof. R. G. Herd . 


- fL S = 








Boy 
MA LNOMJ 


JAVA O038H3^OD "03 
Aq18NW3SSV 


313 Ἴ4ΛΟὈ HO HILIAS 





αμ, 


YILINONVA 





1115 o=- 





( 

SN 

a 
| 
















8313WOWH3HIL 
1103 HLIM 3A1VA AI 
lU L 1398V1 
MM \_ MBSNVHO (OY Ph; Е 
- auod Al8W3SSV SNI 
- => NI1113 
NOIL23S JVOILH3A -I 13 ον INV 


XO8 A18W3SSV ‘HWO4 'L139MVL 





is attached (with the use of clear glyptal) across the slit in the 
valve, which is then raised into a vertical position. In this position 
the system is aligned so that protons passing first through the beam 
defining slit strike the foil at apvroximately the middle of the 

valve slit. 

If stopping power of a gas is renuired, the foil across the 
valve slit is made gas-tight by painting with glyptal around all the 
foil edges. The system is placed on the vacuum system and the valve 
is raised end clamped against the valve seat to gas-tightness. (The 
valve may be manipuleted from the outside by a shaft through a Wilson 
seal.) Gas may then be introduced to the desired pressure, and the 
experiment performed. This technique interposes both a foil and the 
gas in the path of the beam at the same time, so that in determining 
the stopping power of the gas the energy loss in the foil must be 


known and allowed for. 


В. Preliminary Measurements on Foils. 

The foils used were of copper and nickel, of nominal thickness 
between 0.03 mil and 0.1 mil. This amounts to a surface density of 
from 0.6 to 2 me/cré. This is somewhat thicker thon is usually used 
for stopping power experiments; but in selecting such foils the 
greatest importance was placed on consistency in thickness - the 
lack of which has often plagued experimenters who used thinner foils. 
Furthermore certain fixed errors, such as that caused by carbon 
accumulation on the foil, are less percentagewise when thicker 
foils are used. (See Section IV.A, below.) 


The foils were obtained from the Chromium Corporation of 
tn 





America, Vaterbury, Conn., who prepared them by electrolysis.+ 


They were guaranteed by the company to have an average thickness with- 
in 10$ of the nominal thickness, and to have a spot-to-spot varinticn 
in thickness of less than 57. The foils used were judiciously selec- 
ted to obtain the best consistency possible; and & techniaue was 
developed to check thickness variations to within about 1%. 

This technique depends on the range-straggling of alphe particles. 
The apparatus is quite simple and is illustrated in Figure 5. The 
principle of operation is explained by following the procedure out- 
lined below. 

With ihe foil removed, data is taken of counting rate as a 
function of height of the movable platform from the fixed surface. 
This gives a curve similar to curve (a), Figure 6. The steepness of 
rise of the curve is determined by the range-stragrling of the alpha 
particles and is important in determination of the sensitivity of 
measurement of the relative foil thickness. Now if the experiment 
is repeated with the foil interposed over the 1/32" hole, a curve of 
similer shape, but displaced by amount A, is able to be determined. 
This distance A is the air-equivalent of the foil for alpha 
particle stopping. 

If we set the apparatus at héight H, to give an operating point 
on a reasonably linear part of the steep portion of the curve, we may 
then take readings with various portions of the foil positioned over 


the hole. If the counting rate changes by amount Ac, then this 





1 Information on hese foils was kindly provided by Ir. M. M. Stern- 
fels, Chemical Engineer for the Chromium Corp. 
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APPARATUS FOR CHECKING 
THICKNESS VARIATION IN 
THIN FOILS 


THIN WINDOW (^ 1.5MG/CM?) 
GEIGER TUBE , SUCH AS 
TRACERLAB- TGC-2 
-HEIGHT FIXED- 





Ew TO BE STUDIED 
с NE s. 1/32" HOLE 


THIS APPARATUS = 
MAY BE RAISED OR ——THIN DEPOSIT OF 


LOWERED BY TURN- POLONIUM (~10 JC) 
ING SCREWS 


Ш: 
IT: 
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COUNTS / MIN. 


CALIBRATION 
CURVES 


FOR CHECKING VARIATIONS IN 
FOIL THICKNESS 


CURVE (A), WITHOUT FOIL 





AC 


WITH FOIL 


CURVE (B); 





indicated a change in air-equiuwwlent of Ah . The value of Ah/A 
indicates the fractional variation in thickness from the original 
spot to the new spot studied. If a large number of spots on the 

foil are studied thus, we can determine the mean value of H by 
knowing the overall average counting rate. The statistical variation 
in thickness is easily computed as the averages root-mean-sauvare of 
Ah/A . Such a variation corresponds to a counting variation 
(fractinnal) of Ac/C . To translate a counting percentage devia- 
tion to a thickness percentage deviation, we obviously have to mul- 
tiply the former by 


Factor = &.4h _ C . (III-A-1) 
A Ac ~ fA e(slope curve at opert'g Bt.) 


If a finite number of counts are taken at each point, the count- 
ing rate is subject to a statistical variation due to the randomness 
of the radioactive emission process. Such statistical variation must 
be subtracted from the total variation to give the most probable 
variation due to thickness variations in the foil. Furthermore the 
statistical variation from the radioactive emission vrocess 214998 
a practical límit on the accuracy of this technique. 

Table IV (a) gives a set of sample data taken on one of the foils; 
and Table IV (b) summarizes the results of measurement on all the 
foils. 

We conclude from the data that the standard deviation of the 
thickness measured over spots 1/32" in diameter is about the order of 
magnitude of the precision of the measurements. “e estimate that we 
can cite for all the foils a standard deviation of about 1.4% as 
measure of the spot-to-spot variation from the mean. Thus, if we 





Table IV. 


Spot-to-Spot Variations in Foil Thicknesses 
Spot Size : 1/32 " 


Part (a): Sample Data - 0.1 nil Copper Foil. 


Spot No, Counts See. Cts. /min. Ac (Ae)? | 






1 4096 828,5 297 - 1.4 1.96 
2 " 770.3 319 +20.6 424.36 
Š " 800,0 307 + 8.6 73.96 
4 П 831,0 296 - 2.4 5.76 
5 д 844,0 291 - 7.4 54.76 
6 τ 827.0 297 - 1.4 1.96 
7 i 842.0 292 - 6.4 40,96 
8 de 855.0 288 -10.4 108.16 





AV. = 298,4 711.88 


Standard Deviation (small sample theory) = ¡Eh 


= 10.1 , or 3.4% 








Standard Deviation of Radioact. Decay Rate = 1 
4096 
Р == 1.56% 

C 298.4 Cts./min. 

A š 3,0 teens / Conversion Factor = 292.4 κ6 

Slope = 100 Cts./turn ae 

am 0.37 
Part (b): Summary of Results. 
No. of Total Std. Net Std. Conver- Std, 
Foil Spots Std. Dev. of Dev. in sion Dev. of 
Measured Dev. Radioac. Count.Rate Factor Thickn. 
.05mil Ni-1 9 5.9 2.2% 5,47, 294. 2.97 

t 11 2.7 ο. 3.0 . 50 1.» 

" 13 Seo асв 2„°7 «33 1,4 
eOSmil Ni-2 10 de 8 e ec 3.0 „48 1.4 
.OAmil Ni 5 5,5 ca 2.7 „6? 1.8 

" 9 2.9 Qe 1.9 .7ı 1.4 
.O3mil Ni a 260 E 1.4 ooo ls 
1 са 5 1 5 a 0 - small 
el mil Cu 8 3.4 1.6 5.0 37 1.1 


| 
| 
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use average thickness for our computations, while the proton beam 
passes through the foil at a spot selected at random, a probable error 
in stopping power results of about 1% ex&sts because of thickness 
variation alone. 

We ean draw the above conclusion because the 1/32" spot size 
used is about the same as the diameter of the intense "core" of the 
proton beam fram the Van de Graaff generator when it is well focused, 
We cannot as easily make the same conclusions with regard to strag- 
gling measurements, since variations of thickness within a sincle 
1/32" spot will increase the straggline without affecting average 
energy loss. 

It is possible to make a hypothesis that the varietions using a 
1/32" spot are due to more radical variations within the snot of this 
size which are not completely averaged out. This hypothesis is opposed 
to the hypothesis that the variations are broad and slowly chanzing 
in comperison to the 1/32" distance taken as a unit. These were 
tested by making a series of runs in one of the nickel foils with 
successive spots btudied being contiguous and slong a straight line. 
The results were quite consistent with the hypothesis that the varia- 
tions in foil thickness were quite gradual, and highly inconsistent 
with the hypothesis of more rapid, stronger variations. 

This still does not obviate the possibility that very finely 
spaced variations may exist which average out almost comnletely over 
a 1/32" spot. The possible existence of such varistions therefore for 


6 


z -2 
spot checks of size 10 ~ to 10 centimeters must still be admitted. 


Such variations may cause errors in the experimental results on 
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stopping power straggling and are discussed in the section eoneerned 
with such errors (Section IV.D.2). 

The average thicknesses of the foils, measured in units of mng/em?, 
were determined by weighing them and measuring their linear dimensions. 
The weighing was accomplished by microbalance in the Licroanalytical 
Laboratory, Chemistry Depariment, Ohio State University, through 
courtesy of Professor W. M. MacNevin. The measurement of the linear 
dimensions was accomplished under s low nower ee and checked 
by steel rule and reading glass. The foils were rectangular in shape, 
as checked by assuring oneself of the eouality of the diagonal dimen- 
sions, and therefore the area was taken to be simply the vroduct of 


the two linear dimensions. Results of the foil measurements are ~iven 





in fable V. 
Table V. 
Determination of Foil Thicknesses 
Foil Area (en? ) Weight (mg) Thickness 


(mg /em*) 


1.150 £( .45%) 
1.175 £(.65%) 


.05 mil 111-1 8.10 + .01 
.O5 mil Ni-2 3:06 EOL 


0.640 + (.457) 
1.155 + ( . 3⁄7.) 
x) 2.20 E (WS) 


.03 mil Ni 4.75 4: .Ol 
.05 mil Cu 2 8.65 + .01 
mil Cu 16.30 + ,01 


(. 
(. 
¿04 mil Ní 7.23 + .01 (. 0.977 +(.457) 
(. 
(. 
(. 





We may conclude that the probable errors inherent in the use of 
the thicknesses given are the statisticsl sum of the 1% due to non- 
ugifornity and 0.457 due to error in measurement of average thickness. 


This gives a total probable error of about 1.1% in estimating average 


^7 








proton path length through the fcils by the above values of foil 


thickness. 


C. Design of Magnet Current Regulator. 
It was found upon first operating the Van de Graaff generator 


that the current pwovided to the analyzing electromagnet vas not 
sufficiently stable for the close energy resolution desired. It 
was found that for a given setting of the controls, the magnet 
current tended to swing to either side by a variable amount, giving 
a deviation from the mean of as much as 0,1% . 

This amount can easily be translated to energy variation. If 
we accept a rough formula (see Section IITI.F, below): 

B= (45) , (ттт) 
it is easily shown thet a deviation of 0.1% in the current, I, is 
equivalent to a deviation of 0.2% in the beam energy, E. Thus at 
beam voltage of 1 Mev, it was imvossible to prevent the voltnge from 
wandering ss much as 2 kev to either side of the mean. 

It ts true that the operator has manual controls for setting 
the magnet current; but the variation was so erratic and sudden at 
times that the operator could not anticipate and counteract the more 
rapid swings. The disadvantages inherent in bhis variation are manifest: 
(a) They make it difficult to maintain and read accurately the 
mean setting of the electromagnet current, thus preventing accurate 
calibration of the beam energy in terms of magnet current. 

(b) Since beam position at the slit system (see Figure 3), located 
after the magnetic analyzer, is used to ap LEEN accelerator 


woltage (Smith, 1952), magnetic current variations are reflected 
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in beam energy variations. The beam stabilizer, that, is, can never 
prevyert beam energy variations corresponding to magnet current 
variations. 

In terms of the variables we are trying to determine, this var- 
iation would mean an uncertainty in peek position of about 1 kev 
added to all other uncertainties and вобла add greatly to the f»ctors 
eonstituting the Tia discussed in Section II.B above. It was found 
expedient therefore to provide some measure of regulation of the 
nagnet current. 

A brief discussion of how the magnet current is provided would be 
appropriate. <A motob-generator set puts out a d.c., suonly of maximum 
rated values 360 volts and 4.75 amperes, The output voltage is 
controlled by varying the d.c. voltage to a separately excited field 
winding on the generator side. The d.c. voltage for the field is 
obtained by full-wave rectification of 6 variac-controlled elternating 
voltage, obtainéd in turn from a Sola-regulated 115-volt a.c. power 
supply. Thus the magnet current can be varied by changing the 
setting of the variac supplying the field current. This actually 
varies the voltage from the d.c. generator; a fine control of current 
is obtained by means of a variable rheostat in series vith the magnet 
coils. 

A capacitance of 200 Mid. is connected in parallel with the 
magnet. Since the magnet coils have a high inductive reactance 
(inductance estimated to be about 30 henries), fluctuations of 
frequency much faster than one cycle per second pass through the 


condensers rather than the coil. Aside from this, and the Sola 
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regulation mentioned above, no other regulation wes in the circuit 
prior to the time the regulator herein described was installed. 

Several possible schemes for regulation of current itself were 
considered and rejected. It was decided that uncontrollable resis- 
tance changes in the magnet circuit (caused by temperature changes in 
the coils) were slow enough to be compensated for manually by the 
Operator. Furthermore fine control of current by a variable resistor 
was desired. ‘Thus, a voltage regulator was actually used for stebil- 
izing the current. The d.c. generator voltage is affected not only 
by the field current, but also by the motor speed, which may nossibly 
experience small variations. It wes decided therefore to put the regu- 
lator directly in the magnet circuit between the generator andá the 
load, which consists of magnet coils and the control rheostats. 

It was necessary under the cireumstances to keep the svstem 
quite simple, since any complex electronic apperatus has a very low 
power efficiency and the output of the d.c. generator is rather limited, 
The very simplest type of regulator of the "transconductance" type is 


used (Hunt and Hickman, 1939). 


+ + 
| Ry | 
E ra 
| | DC 
| 


4|'- | Output Rr, 


-- — — p — = А 


Figure 7 shows in shhematic fashion the regulator in »rinciple. 
The triode is in parallel with the load, and the current through ít 


adds to the magnet current to give the total current required from 
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the d.c. source, The amount of tube current depends on the output 
voltage and the grid voltage. If a small positive variation occurs 

in the input voltage, more current flows fram the source. This alee 
also causes a rise in grid voltage, which produces a greater current 
drain through the tube. “ith ppoper ehoice of the circuit parameters, 
the current increase through the tube matches the totel current rise 
from the source - causing the same amount of current through the loed. 
Or, from another point of view, a current rise through the tube 
creates additional voltage drop across Ry |, which is designed to 
just counteract the original voltege rise at the input side - this 
giving the same output voltage. 

A rigorous proof of the above statements, indicating how the 
circuit psrameters are selected is given in the anvendix. Also in 
the apvendix is a circuit diagram of the complete magnet current 
control snd regulating system, as it exists at the present time, with 
further explanatory comments on its operaticn. 

Results with the regulator have been excellent. With prover 
attantion on the part of the operator to compensate for slow 
resigtance changes and keep the regulator controls properly adjusted, 
the current variations are cut down to sbout one-tenth their value 
without regulation. This means a voltage varistion of the beem of 


about 100 - 200 ev , which is small compared to other variations and 


errors. 


D» Technique for Taking ana Correcting Yield-~Energy Data. 
With the chosen target in place and a vacuum on the whole svstem, 


the Van de Graaff generator is started and protons produced in 
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accordance with standard procedure. The magnet is then "cycled" 
twice by running the magnet current up to a relatively high velue 
(about 3 amperes) and then back to zero. ‘The subsequent data is 
taken with ever increasing magnet current, so that one stays consis- 
tently on a reproducible magnetization curve (magnetic flux versus 
magnet current), without any hysteresis effects. 

A8 and 1f required, the stopping foil is inserted and clamped 
at this point in the ppocedure, and the desired amount of gas (which 
may be zero) is introduced into the holder (Figure 4.) . 

With the magnet current set at the initial (lowest) value, the 
accelerator voltage is set so as to send the protons through the 
defining slits. The proton accumulated charge is counted by an 
"integrator" with a meehenical register and the gamma-rays detected 
by the usual Geiger counter-scaler apvaratus. At each current set- 
ting the data taken includes: magnet current, total gamma-ray counts 
(usually set to a predetermined figure), time to collect these 
counts, total integrator count of accumulated proton charge. Arrange- 
ments are such that the proton integrator operates only while the 
Geiger counter-scaler is accepting and registering gamma-ray counts. 

The gamma-ray counts must be corrected for "background", This 
background is a function of the generator-accelerator voltage, since 
x-rays from the generator provide a sizable proportion of said back- 
ground. Experiments mode without a target indicate that accumulnted 
background can be estimated by dividing the time (seconds) for the 
prescribed number of total counts in e single run at a current set- 


ting by a factor, given in Table VI. These factors are valid only 


— 
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Table VI. 
Factors for Determination of Packground 
When Operating Van de Graaff Generator 
Type Tube: Technical Associates Beta Counting Tube TA-Bl. 


Location: As close to the target as possible - separated from 
sane by thickness of target holder, 1/8" of lead, and 
about three millimeters of air. 


Position: (One tube) - Side to target, with center line at mid- 
target height. (Two tubes in parallel) - One above the 
other, sides to target, with mean of positions of 
centerlines at mid-target height. , 

Shielding: Tube(s) mounted within a two-inch thick lead box, with 
hole in same just large enough to permit ready inser- 
tion of target holder. 


ee ee ee ee eo «αρ 





Analyzing Magnet Factor by !ihích to Divide Length of 
Current (amp.) Run in Seconds to Get Background Counts 
One Tube Two Tubes in Parallel 

295 3.0 2.0 

. 60 2,8 2.0 

$50 2.6 2.0 

. 70 2.4 2.0 

220 2.2 1.9 

„80 260 кү 

„85 1.8 1.5 

. 90 1.9 1.4 

„95 1,2 p 
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1f the geometry, tube type, and other factors are reproducible; and 
these variables are indicated in the table. 

The proton integrator count for a single run is obtained by sub- 
tracting the register values before and after the run. The values thus 
obtained also require correcting, as a certain amount of charge leakage 
oceurs through paths other than that which actuates the register. By 
experimentally checking the integrator register counts ver unit time 
against galvanoneter readings of current fron the Faraday cage to the 
integrator, it hos been found that integrator readings must be increased 
by an amount corresponding to about 35 eounts per minute of time for 
the run, Such counts, when corrected, thus provide a more constant 
ratio between the figure for total integrator counts and the actual 
number of protons striking the target. 

The "yield", on an arbitrary scale, is the result of dividing 
the corrected garma-ray counts by the corrected proton integrator 
count. 

These data are taken at successively higher msgnet current 
settings. The desirable current increments between successive 
readings for our experimental work was as follows: for sharp 
veaks, current increments of 0.5 ma. are necessary; for broadsr 
peaks, data every 1 ma. are sufficient. 

If gas stopping is being studied during any senuence of indi- 
vidual runs, the manometer, thermometer (see Figure 4), „nd time of 
day are recorded just before ges introduction, just after gas intro- 
duction, just before the gas is finally removed, and just after the 


gas is removed. Readings at intervals between start and finish may 
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also be made if desired. 

The yield can be plotted against magnet current, or if the 
analyzing current-proton energy relationship has been established 
(called "calibration" herein) the yield can be plotted versus nroton 
energy to give a more proner spectrum. For our purposes only the 
first type of vlot has been recuired, as it is possible to read the 
veak maximm values and the statistical measurements of width 
directly in terms of magnet current and then convert t° the proton 
energy scale. 

Table VII gives a sample set of data for a single peak, with 
yield computations, Figure 8(a) shows the plot of the resonance peak 
(unshifted) for this particular set of data; Figure 8(b) is a plot of 
the same peak, shifted by insertion of stopping material. The 
probeble error of the vosition of each voint is indicated by the 
size of the crosses. The error in horizontal direction is based on 
personal observation and judgment; the error in the vertical direction 
is determined from strndard statisticz1 deviations of the total num- 
ber of gam-ray counts (P.E. = 0.67 VNo. of Counts ). In the latter 
case, the actual probable error is somewhat greater, but the largest 
errors come from counting statistical errors and give © rough basis 
for estimating the extent to which our smooth curve has to follov 


the experimental points. 


E. Target Preparation. 


The target noterials used were either lithium fluoride or 
aluminum. The targets were formed by coating a blank tantalum 


disc with a thin coating of the target substance. ‘The coating was 
- Z4 


Table VII . 


Sample Set of Data for Computation of Yield and Plotting of Peak 


Peak Covered: 340.4 kev - FL9 , Target Used: LiF-2 
No Stopping Medium Interposed. 


One Geiger Tube Used, 






— -- 


Magnet |Proton | Time | Correc. | Net 
Current| Int. (see) 





(amp) ICounts Counts |Counts 
.4900 49 28.6 17 66 
10 82 45.0 26 108 
20 83 56.6 35 116 4 
50 59 34.0 20 79 64 55 „7? 
40 94 57.2 33 127 128 18 110 „9 
50 84 46.8 27 111 128 15 115 1.0 
60 80 45 ,4 25 105 256 14 242 2.) 
70 69 41.8 24 93 256 13 CAD 2.6 
80 61 37.7 oe 83 512 12 500 6.0 
90 62 44.0 26 88 1024 14 1010 LISS 
~ 0000 43 32.0 19 62 1024 10 1014 
10 71 38.4 22 93 . 2048 12 2036 
20 64 37.4 22 86 2048 12 2036 
30 105 60,0 35 140 4096 19 4077 
40 136 του 42 178 4096 23 4073 
50 57 33.6 20 77 2048 1] 2057 
60 55 95,5 2] 74 2048 11 2037 
70 83 53.6 31 114 2048 17 2031 
80 61 38.0 22 83 1024 12 1012 
90 70 46.2 27 97 1024 14 1010 
«9100 83 62,6 37 120 1024 20 1004 
10 224 84.0 49 273 512 26 486 
20 85 95.0 32 117 256 17 269 
30 67 39,0 23 90 128 12 116 
40 77 45.6 ZO 102 128 14 114 
60 70 41.5 24 94 128 15 115 
80 110 6237 37 147 128 20 108 


These results are nlotted in Figure 8(a). 
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Table VITI. 
List of Resonance #nergy Values for Proton-Gamma Process in 


Materials Used as Targets 









Element Reference Resonance Width at 
Energy (kev) | Half Ht. = Í sa 
(kev) 





Li” (1) 441.1*.5 12.2+.5 
Е19 (1) 540,4 2,5 
598 37 
669 7.9 
OTOIO ο ο 
985.9 8,0 
p" (2) 505 From 1 
630 ev at 
771 600 kev, 
986 to 1 kev 
1112 at 1.4 
EOS ev 
1176 
References: (1) Ajzenberg and Lauritsen, 1952. 


(2) Brostrom, Huus, end Tangen, 1947, 
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performed by a rather standard technique of evavoration in ϱ 
vacuum (Strong, 1945). The materials selected were those for 
which the energy resonance levels and peak widths sre known with 
reasonable асеюгасу. The levels for lithium, fluorine, and 
aluminum as obtained from the latest and best references are 
tabulated in Table VIII. 

A succession of targets were made and used, before sone were 
obtained which were entirely satisfactory. A satisfactory target 
must have the following qualities: 

(a) The target should be thick enough to give high yields of 
gamma-rays per proton count, thus providing good counting "statistics". 
(b) On the other hand, the target thickness should be small enough 
that target thickness, T, ana fl, (Section II.B) can be minimized, 
thus minimizing the effect of errors in their estimetion. nd under 
any cirewnstances T must be much less then the energy difference 
between adjacent peaks, to avoid interference fron other peaks when 

@ certain one is analyzed. 

(c) A fair uniformity of target thickness should be provided, so 

that when a proton beam is somewhat diffused over the target, the 


same values of T ana f), will apply to all portions of the bean. 


F. Analysis of Auxiliary Statistical Variables. 

Preliminary to the primsry calculetions of stopping power sre 
those calculations of the other variebles which are a part of edun- 
tions (Ii-B-4) to (II-B-7), inclusive. These calculations require 


a knowledge of the calibration relation between analyzing magnet 
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eurrent and proton energy. Unfortunately an accurate determination 
of this relation requires in turn a knowledge of target thickness, T, 
which is one of the variables referped to above. Originally a rough 
assumption as to the proper calibration relation had to be assumed 
and calculations made by the process of successive approximations, 

In this report we may start with reasonably accurate information on 
the mean calibration curve, which has actually been obtained at a 
later stage in the experiment (See Section III.H). 

Within a limited range it is found by plotting calibration data 
on log-log graph paper that the calibration relation can be expressed 
by a formula of the type: 

Bes kr, (ral) 
where E represents proton beam energy (Nev) anà I is analyzing 
magnet current (amperes). (See Figure 9.) 


Besides K and a, we wish to find the value of F, where 


эш dE - al 
T > dT = ак 


which is used to translate magnet current variations into energy 


А (III-F-2) 
variations at some particular peak. It has been found, using data 

to be presented later (Section III.H), that the values for K, a , and 
F given in Table IX are fairly accurate. 

Let us now analyze sme of the statistical distributions in 
energy a little more closely. It is necessary for our purposes that 
the ststistical spread in these distributions be expressed in terms 
of the standard deviation about the mean or its square, the variance, 
Many of the distributions are of such a type that the spread is not 


usually expressed thus, and calculations must be made to permit 
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Table IX. 


Values of constants К, а, арпа F (<= dE/dI) in the 


Magnet Current-Proton Energy Relationship: E καπ. 
Current Approx. 
(amp) Energy > E > 
(Μεν) 
¿50 344 1.355 1.975 1.362 
250 „417 1.353 1.972 1.490 
«60 „494 1.550 1.970 1.620 
e GS .579 1.348 1.967 1.748 
e 70 „669 1.545 1.960 1.870 
.75 ‚765 1.339 1,950 1.986 
«80 .868 1.338 1.935 2.09 
.85 e 972 1.999 1.920 2,19 
e 90 1.076 1.335 1.905 PB 
„ 95 ΕΤ 1.505 1.887 3255 


translation to those statistical variables we employ іп common for 
all the distributions. Resonance curves for example have their 
spread expressed in terms of a width of peak at half the maximum 
height. 


1. Method of Obtaining Sl... 





As just noted, it is necessary to translate the natural peak's 
width at half height, designated as Ps ¿ » to the standard deviation 


101 


since the theoretical computation of Mi takes one into compli- 


nat’ 9 8τ6 immediately faced with a somewhat difficult situstion, 
cations of both a mathematical and a practical nature, 

The shape of the theoretical resonance curve is given by the 
well known Breit-‘Jigner formula (Breit and Wigner, 1936), modified 
for the proton-gamma reaction (Bethe, 1937): 


m. 
mE. —— 








an | u , (Ттт-Р-з) 
ЇЕ (E - Ees) +f Ko 


where O is the probability cross-section for an atom to undergo 


Oo = 


the proton-gamma reaction. oe is the barrier penetration fector 
for the proton in trying to pierce the Coulomb potential barrier of 
the atom, and at high proton energies (compared with Coulomb barrier 
height) it is dlose to unity, whereas at low proton energies it 
equals approximately exp (-K'/VE ) š 

This expression is difficult to deal with mathematically, but 
we may make further simnlifications for the sake of discussion. Since 
most of the resonances used are extremely narrow compared to the proton 
energy (E) exciting them, the factor aS is very slowly varying in 
the immediate vicinity of the peak, and can be taken as constant. ve 


then have a resonance peak with a shape given by the formula: 








с х К . (III-F-4) 
2 
(E - Een)! f K, 
р 
Ву definition, Sl, is then given by 
Er 
K dE 
2 
+ = 
1 Ko/(E E ορ) 
2 E. 
ESE š (III-F-5) 
nat E. 
K dE 
2 
(E = Ea! t K 
Ez 


Tbe difficulty arises vhen we try to deside vhat limits to take for 


E, and E, , on either side of the peak. E_ cannot go below gero, 


but there is no natural upper limit to select. The value of the 


a 
P. 


cu 





denominator is finite with limits zero to infinity; but the value 
of the numerator increases without limit as E, aprronches infinity. 
Thus the value of the standard deviation will depend upon where we 
choose to take this upper limit end befoames somewhat arbitrary. The 
inclusion of the factor c ҮІ does not simplify our task, since 
the numerator still diverges with increasing value of ER. 

Et 

To avoid the decision as to where precisely ve cut off the 
"tails" of the distribution, we shall employ the same apnroxímation 
used by Madsen and Venkateswerlu (1948b). Noting thet the resonence 
curves are rather similar in shape to a Gaussian, or "normal", dis- 
tribution curve, one may assume that the ratio of the standard 
deviation to the width at half-msximum height is the same for a 
theoretical resonance curve as for the true Gaussian shaped curve. 
From tables of the latter function it is easily found thet as a 
consequence 


n OE ы - (III-F-6) 


nat ^ 
Table VIII, vhich gives the resonance vidths for the target 
materials used, also gives the values of (at conputed by the above 


formula. 


2. Method of Determining (1, В 


The distribution curve for energy loss in the target materíal 
is quite simple. Ye may assume that the roton stopping power is 
practically constant over the small energy variation the beam under- 
goes in the thin target, so that this distribution curve has the 
shape of a rectangle. That is, if y( AE) is the frequency variable, 


then it is defined as: 
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0, for AE,<0 am Åp >ti 


P" = > (III-F-7) 
k, for ο < AE, = Т ° 
Then, À : T 
D - кт k (OE = м)? d(4E) , 
0 
m ure, (III-F-8) 
so that 
N, = „288 T : (III-F-9) 


3. Determination of N, б 


The distribution of energies of protons passing through the 
slit cannot be known with the degree of certainty one would like. 
Experience indicates that for a narrow slit the beam energy stabili- 
zer can usually hold the center of the beam within the edges of the 
slit quite well. The maximum variation in proton energy which will 
still permit the beam center to pass through the slit has been can 


puted by Smith (1952) to be: 


ome > AR : (III-F-10) 


where w is the width of the slit in millimeters, and ER. end E 

are in the same units. This formulae is probably not highly accurate 

in view of the assumptions and apvroximations involved in the derivation. 
To provide a check, the shift in location of a peak at 669 kev was 

determined experimentally when the slit position was shifted by 4 

millimeter. The shift in peak position was 1.2 ma, of magnet current, 

which corresponds according to Table IX to a beam energy shift of 


1.87 x 1.2 = 2.25 kev. These values are consistent with an equations 


A ae 





Е м 
dis m 150 ° (III-F-11) 





This experiment probably does not give a highly accurate 
answer, but it is considered an improvement numerically over the 
theoretical formila's prediction; and we shall use equation 
(III-F=11) therefore, 

If we assume that energy variations in the accelerator are 
such that a proton is equally likely to pass through the slit with 
any value in the corresponding spread of § ES » we have a rectangu- 
lar distribution; and as in equation (III-F-9) we obtain: 

M - .288 Ew/ 150 = Ем/ 520 . (III-F-12) 


Since for all the experinents herein described, v equals one 


millimeter, the formula becomes: 
(L. = 5/520 , (III-F-1>) 
where the units of E and fis are the same. 


4. Determination of fly. 





The factor n. is impossible to obtain directly, and it can 
only be guessed at vhen the resonence peaks vidths ere anoelyzed, 
There is no reason to believe that it is even constant from one 
day to another under what are presumed to be similar operating 
conditions. Since current variations cause greater beom voltage 
variations with higher energies, we may guess that f), will general- 
ly increase with increasing proton energy, though this is by no 


means certain. 





G Analysis of Unshifted Resonance Peaks and Determination of 
Target Thicknesses, 

Table X gives the analysis of the statistical variations of 
most of the resonance peal:s studicd without stopping nedia interposed. 
This analysis is based upon formula (II-B-5), and shows all camputed 

2 2 
values Tor ir . Observation of these results leads to the 
following conclusions: 
(a) There are connaratively large uneertainties present in all of 
the data presented, so that the final probable error in the estina- 
tion of In... is on the order of 1 kev or mere, This is especially 
evident from the size and number of negative values obtained, since 
all values shovld be positive for variances. 
(b) There seens to be little reason to believe that Lt, is large 
enouch to consider, except in a few cases. We shall. therefore imore 
it in our computations but must consider its nossible influence when 
the matter of errors im results is discussed. 
(ο) For low energy values N is at a maximm, vhereas SL, is at 
a minimum. Therefore in determination of target thickness, greater 
weight should be siven to data obtained from lower energy peaks. 

In ofder to exoress the target thickness in terns of stopping 
power, we rust standardize sch values to a specific value of proton 
energy (1 !iev is usıal). This means that after having determined 
target thickness at a certain proton energy we must convert this 
value to thickness for a 1 Mev proton beam. This сап be done if we 
know the relative stopping pover of the tarzet substance at І liev 


compared to the stopping power of the target at the energy of the 
- HZ 





Table X. (Two Pages) 


Analysis of Contributions to Variances in Unshifted 


Resonance Peaks. 










8-29-52 Lir-1 540 4 «50 1.562 
y И 441.1 „566 1.553 
п ч 669 «70 1.870 














9-24-52 AL-3 503 „60 1.620 
" m 650 „68 1.820 
" " 986 .85 2.170 
" " 1112 „91 2.255 


9-26-52 A1-5 650 . 68 1.820 


10-16-52 Lir-1 669 . 70 1.870 
H s 873.9 «80 2.086 


" " 935,3 .83 2.1955 


10-17-52 А1-5 650 . 68 1.820 
τ " 986 „86 2.185 
M = 1112 sa сое 


10-25-52 LiP-1 441.1 e066 1.533 


12-2-52 L1F-2 340.4 «50 1.362 
H " 669 «70 1.870 
H H 875.5 „81 „105 
ч H 9355.5 ο84 2.155 


12-18-52 LiF—4 441.1 .57 1.540 
M " 669 .70 1.870 
ч M 873.5 «805 2.095 


2.144 


- £3» 


Table X. (Concl.) 























—— g 

(kev?) (kev?) (kev^) 

Table VIII | (III-F-15) (ІІ-Р-5) 
1 1,90 2,60 6.75 0,94 0.43 9.38 
2 4.95 7,60 57.7 26.0 0.72 31.0 
3 35.94 6,25 59.1 949 1.66 27.5 
% 0.80 1.30 1.69 Negl. 0.94 0.75 
5 0.90 1,64 6.69 Ç 1,47 1,22 
6 0,84 1.82 3.51 0.01 3.60 0.30 
7 1.35 3,05 9.30 0.01 4.58 4.71 
8 1.1 2.00 4.00 Negl. 1.47 2.09 
Э Hop 5.10 9.61 E 1.66 - 1.95 
10 L 2,50 6.25 4,84 2.82 = 1.41 
3.24 1.95 
2, 94 
0.51 
0.79 
1.40 

22,4 

16.2 

51.5 

εν 

- 0,24 

= 4.43 

- 1.36 

Ten 
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Table XI. 


Computations of Target Thickness (kev) 








ps 
* 

e 
“J 


2003 «710 0.75 0.866 3.00 

.630 . 794 1.28 1.105 3.85 

. 630 «794. 2.99 1.59 5.51 

.630 „794 2.94 1,714 5,95 

. 986 „99% пера, = = 

. 986 . 993 0.51 0.714 2.47 
1.112 1.055 4.71 2.397 7.53 
Г 12 1.055 0,79 0.889 3.08 


=F js 





particular peak analyzed. 

Bohr (1948) indicates that within this energy region which ve 
are concerned with, the stopping power is inversely proportional to 
the square root of proton energy. To the extent that data is 
available in this region, such dependence is roughly confirmed 
(Madsen and Venkateswarlu, 1948a; Kahn, 1953; and Figure 11 herein). 
In view of the approximate nature of the data given in Table X, it 
is sufficiently accurate to convert target thickness to the 1 Mev 
standard energy by multiplying the target thickness by the square 
root of the proton energy at which the specific measurement was made. 

Table XI gives the target thickness analysis, with the aid of 
the data derived in Table X. In view of the approximate nature of 
the final results, it is desirable to use other means to aid in the 
selection of definite — If the targets are thin the naximum 
value of the yield at e certain resonance peak will be a rough 
measure of the relative thickness of the target. Thus a comparison 
of the maxima of the various peaks for the various lithium fluoride 
targets will add additional information on which to base final judg- 
ment. It is expecially desirable to do this for determination of the 
thickness of target LiF~4, since it is evidently so thin as to make 
detection from an experinental peak width practically impossible. 
Table XII gives this data, as taken from experimental results, The 
table indicates that in a rough sort of way, the relative yield 
fram the various lithium fluoride targets are in the ratio: 


= 7 / E. 





Table XII, 


Resonance Peak Maxima for Lithium Fluoride Targets 


Maxinum Yields 


340,4, Кет 
441.) 


441.1 
669 


669 
# ©2735 
#o35,3 





* Evidently the target Lif-1 has a somevhat variable 
thickness. We can speak only of its averaze value. 


it Counting zeonetry somewhat variable. 


Table XIII, 


Target Thiclnesses (kev) for Various Proton Enersies 


= 
e 
e 
O 


NOW WW m „ә 
9 


9 
ο ο ὃ- ο δν» 


О 
1.8 
126 
1.2 
132 
1.0 
1,0 





On the basis of all information available, we can now make an 
estimate of target thicknesses for all targets, Table ZIII gives 


- — 





the thickness estimated at 1 Mev, and also gives the computed values 
of target thickness ai various resonant energies. In the computation 
the target thickness 35 assumed, ab before, to be inversely propor- 


tional to tho square root of the energy. 


H,  Encrgy-Current Calibration, 


Since we can now find the value of ES for each peak measured, 
by use of equation (II-B-4), we may find the calibration curves for 
accelerator beam energy (Е) versus analyzing magnet current (I). 
All data from experiment and previous analysis is listed in Table 
XIV in such a fashion that the calibration curves can be put on a 
logarithmic scale, The curves themselves are depicted in Figure 9. 

Table XV is a list of data similar to Table XIV, except that 
the data are for single resonance peaks and are not taken as part of 
a complete calibration series. Their use is indicated below. 

In usin; and interpreting the calibration curves, several 
significant points have been discovered: 

(a) If the data is plotted to a logarithmic scale, vertically and 
horizontally, the calibration line draw is almost a straight line. 
this indicates that an empirical formula of the type E = K Í 

is approxinately valid over the whole energy range, and is quite 
accurate within smaller limits. (This justifies equation (III-F-1).) 
This is a change from the previo’s custon (Grove, 1950; Taylor, 1952; 
Russell, 1952), by which data was plotted as JE versus I and the 
best straight line drawn, This older type of plot presumes that an 
empirical formla of the type E = K (a. is valid. Plotting 
of a typical set of calibration data from Table XIV according to the 
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Date; 
Target; 
[511% 
Setting, 
mm. 


8-29-52 
LiF-1 
4,25 


9-26-52 
А 1-5 
4.75 


10-16-52 
LiF-1 
4.75 


10-17-52 
A1-3 
4.75 


12-2-52 
LiF-2 
3.75 


12-3-52 
A1-3 
3.75 


Enat 
(Mev) 


05404 


08735 


Table XIV. (Two Pages) 


Calibration Data 
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r 


(lev) 


. 0034 
0030 
0024 
«0021 


0021 
0019 
«0015 
.0014 


. 0024 
¿UNE 
.0021 


¿0019 
. 0015 
.0014 


0105 
0073 
» 0064 
0062 


0021 
„ООЗ 
«0015 


— p o win аР --- l —  — wA )À— Ue — - - 


Eo 


(Liew) 


‚5458 


«6714 
.3756 


2051 
‚6519 
. 1875 
1.1134 


«6714 
.8756 
. 9374 


«5519 
9875 
1.1134 


05007 
„6765 
‚8799 


2091 
| 0519 
ου 


1οσ Ε 





| elel ei 
Ф 9 
Sg 

de 

-ᾱ ΟΥ 

‹л © 


034 
006 
945 


η 


9 


4 
O ιοςο- 


I log I 
(amp) 

4975 1,6966 
A 
‚6996 1.8448 
‚8020 1.9042 
„6052 1,7805 
731 1.8201 
.8563 1.9326 
29085" 1.9558 
‚Ross 1.8465 
‚8030 1,9047 
„8315 1,9199 
‚6789 1.8318 
8575 1.9531 
105 1.9592 
‚5036 1.7021 
72036 1.8475 
«80734 1.9071 
.8377 1,9231 
“ο 1.7864 
.6858 1.8549 
‚8640 1.9365 


“Qo Qm ῳ- amas qg Se a сїй Ч= сю Qp. ee eee eee EO ee ӘС ЧӘ eee s Wr 


Slit setting eruals reading of nosition rierometer minus one- 
half the rerding of the slit widtr ricroneter. 


It essentially 


inäicstes the vosition cf the midele of the slit. 








Table XIV. (Concl.) 
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Date; 
Target; log Fo I 
Slit 
Setting, (amp) 
mie 
12-18-55 4411 «0008 ¿4419 1,6453. .5690 
LiF-4 „669 «0006 .6696 1.8258 .6995 
ΞΕ .8735 .0005 .8740 1,9415 „8052 
.9353 «0005 .9358 1,9712 . 8334 
12-19-52 .503 „0021 .5051 1.7084 . 6095 
41-3 630 .0019 «6519 1.8006 „6817 
4„25 22901 ορ .7727 1.8880 . 7562 
«986 „0015 „9275 1.9945 „8620 
1.112 .0014 1.115%" 0.0067 „9176 
1-9-55 4411 «0008 «4419 1.6453 „5700 
Lir-4 .669 «0006 .6696 1.8258 .7028 
3.50 «9705 .0005 .8740 1.9415 «5074 
.9353 «0005 «0555 1.9712 „8567 
1-19-55 «669 «0006 .6696 1.8258 .7014 
LiF-4 «8755 «0005 „8940 1,9415 „8080 
ποῦ .2353  .0005 69358 1.9712 „8567 
1-22-53 (an) .6559 «0006 „6696 1.8258 „7040 
LiF-4 «87965 «0008 8740 1.9415 „8072. 
5„50 
1-22-55( pn) „669 . 0006 „6696 1.2258 „7014 
Lir-4 «8755 «0008 .8740 1.9415 .8056 
5,50 «9555 .0005 .9358 1.9712 .8335 





Table XV. (Two Pages) 


Data for Single Peaks (Unshifted) 


Date; 
Target; 
Slit 
Setting, (Mev) 

mm. 





10-23-52 > 
Lir-1 «4411 0030 „4441 1.6475 .5670 1.7536 
4.75 


10-28-52 _ P 
LiF-1 . 9353 .0021 .9374 1.9719 .8320 1.9201 
4.75 


10-31-52 _ _ 
LiF-1 9355 „0021 .9374 1.9719 8303 1.9192 
4.75 


1-12-53 (am) a = 
LiF-4 „669 „0006 „6696 1295858 . 7009 1.8457 
3.50 


1-12-53 (m) А = 
LíF-4 .669 . 0006 .6696 1.8258 . 7001 1,8452 
3200 


1-13-53 .669 .0006 .6696 1.8258 .6998 1.8450 
LiF-4 
3.50 


1-14-55 _ И 
LiF-4 „669 .0006 „6696 1.8258 .7002 1.8452 
3.50 


1-15-53 Е B 
LiF-4 „669 «0006 .6696 1.8258 .7023 1.8465 
3.50 


1-19-53 τη = 
LiF-4 „669 . 0006 „6696 1.8258 „7024 1.8466 
3.90 


in ee ee Si мз m= 





Date; 
Target; 
Slit 
Setting, 
mm. 


Table XV. (Concl.) 














Figura 9. 
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previous custom and comparing such a plot with the double~losarith- 
mic plot has indicated that the latter plot has less curvature than 
the corresponding plot by the older method. In view of the fact, 
however, that a slight curvature does exist for the calibration 
curves, even on the logarithmic plots, it was considered more 
accurate to use the curves directly rather than to use a straight 
line formula giving the closest fit over the whole energy range. 

The empirical formula (IiI-F-1) and Table IX are to be used only for 
computations at localized energy regions. 

(b) Different calibration series give curves which are generally 
slightly displaced from one another, although they are essentially 
parallel. This indicates that personal errors during a series of 
runs are small compared to those discrenancies between 

different daily series. Every effort was made in doing the experi- 
mental work to use a standardized technique in oneration of the 
accelerator and auxiliary apparatus; but some variation from day 

to day seemed unavoidable. Such variations nay be due to: (1) lack 
of precise control as to the horizontal and vertical position of the 
beam axis (as indicated by the neutral beam, undeviated by the 
magnetic field); (2) lack of control over the temperature of the 
namet itself. The first effect means that inhomogeneities in the 
magnetic field may be sufficient to be noticeable. The latter effect 
is related to the fact that magnetic elements in an electromapnet 
are more easily aligned when the iron is at a higher temperature, 

so that the same current gives a slightly greater deviating flux, 


The fact that different curves are closely parallel means that once 


PA < 
— ΦΝ 





a few typical calibration curves have been plotted to sive the 
general trend, further curves may be indicated by having data for 
only a single calibration point on the curve. This fact has been 
utilized in the latter part of the research, and Table XV vives 

the data for these single calibration points. 

(c) It is of some interest to note the actual variation of the 
calibration data for one specific valve of energy, to give an idea 
as to the spread in the calibration curves over a six months! period. 
A check of Tables XIV and XV shows that for the 669 kev fluorine 
resonance peak, by use of the same target and the same slit position, 
the average value of the analyzing current was .701/, amperes, with 
extremes of „6696 and ,.7040 amperes. The standard deviation of all 
readings from the mean is .0014 amperes. If on the other hand, one 
were given the current at about this point, the extreme energy variation 
corresponding to this calibration curve variation would be 8.2 kev 
(determined by using factor F, Table IX). The standard deviation 
in energy comes out to be 2.6 kev fro the mean calibration line, or 
a probable error, if one used the mean line, of about 1.7 kev, This 
fact indicates that if ve have some additional information as to the 
proper position of the true calibration curve, sich as by a single 
calibration point, a much better choice can be made for energy 
correspondin;r to a certain value of magnet current; and the probable 
error in such a case would undoubtedly be substantially less than 1 
kev. It should be noted that this conclusion is justified only if 
tho technique is used for data taken on the sane nammet cycle as the 


calibration point. A recycling of the magnet may require a new curve. 
(For exarmles, sec Table XIV, data for 1-22-53, a.m. and p.n.) 
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IV. Experinental Results on Stopping Power. 


A, Possible Sources of Systematic Error and Corrections Required. 


When the experimental data on stopping are analyzed, the possi- 
bility of certain consistent factors requiring correction mst be 
considered. | 
1. Nuclear Collision Losses. 

It can be visualized that proton-nuclear collisions, though 
elastic, will result in a certain small amount of energy loss by 
the proton in the laboratory system. The usual formulas for 
stopping pover do not include this effect, and therefore in theory 
the nuclear sto ping effect should be subtracted from experimental 
data in order to give a result which can be compared with theory. 

Bohr (1948) has derived a formula for nuclear stopping, given 


as follows: 


QE = Wax ame j 2 (IV-A-1) 
Mm v^ G 
2 
where б 2 27 24/3 m [Vo : (IV-A-2) 
М у e 


In order to determine the magnitude of this effect, an example 
was calculated for hypothetical circumstances conducive to a large 
answer. For energy loss due to nuclear collisions of protons of 
mean velocity 0,8 x 10” cm/sec (energy of 332 kev) passing 
throuzh a 0.05 mil nickel foil, the calculated value of A, E is 
0.171 kev. It can be secn by comparison with electronic stopping 


enerzy losses for the foils (Figure 10) that the nuclear effect is 
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only on the order of 0.1% of the total stopping effect. The nuclear 
stopping may therefore be neglected. 
2. Increase of Path Length by Scattering. 

Another effect, also due primarily to nuclear collisions, is the 
deviation from the beam axis of many protons, especially for the case 
f slow protons through large thicknesses of matter, This effect is 
observable visually during stopping experiments with gases, The ex- 

citation of the gas atoms by the proton beam brings about the emission 
of light which marks the proton paths. It can be seen that the bean 
does indeed diverze in a very narrow cone from the entrance foil to 
the gas chamber. Such effect increases the average amount of matter 
traversed by the beam. The geometry of our chamber (Fisure 4) is 
such that the target subtends at the entrance foil an angle of 
5.6 to each side of the axis. Any proton deviated so much as to 
pass outside this cme does not hit the target and cannot contribute 
to the gamma ray yield and thus has no effect on experimental values 
of stopping power. Only at very low proton energies (as noted visu- 
ally) does there exist an appreciable deviation of the protons away 
from the axis. Even in such circumstances the correction for longer 
path length is only a fraction of соз 5, 69 › or less than one-half 
of one percent, At higher energies the correction is still less. 
Since this correction is much less than the over-all experimental 
error expected, ve shall ignore it. 
3, Carbonization d Foils. 

There is a tendency for all parts exposed to the beam to accu- 


mulate a dark deposit of what is called "carbon", though it 
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probably includes a moderate percentarze of hydrogen, This deposit 

is a well-lmown problem for many accelerators. It is presumably a 
result of deposition of volatile hydrocarbons on surfaces plus sub- 
sequent molecular breakdown by proton bombardment and by localized 
heating. It is even possible that if a fol or tarzet is heated to 

a high temperature, hydrocarbon molecules in gaseous form may be 
decomposed upon random collisinynwith the spot heated by the bean, 
even if there is no tendency otherwise to adhere as a semi-permanent 
or permanent film. The sources of the volatile hydrocarbons are seve- 
ral, the nost probable ones being: 

(a) vapors from the oil diffusion pump of the vacuum systen; 

(b) volatile mods from organic substances к as vacuum rings and 
gaskets in the systen; 

(c) vaporized solvents from freshly applied sealing paint (such as 
glyptal) used to stop small leaks; 

(d) vapors from small amounts of acetone or similar penetrating 
fluids used for leak testing; 

(ο) vapors from residual anounts of organic solvents used for 
cleaning parts or the cysten, 

These factors are listed in such an order that the first ones are 
most difficult to avoid, and the latter ones the easiest. Fortu- 
nately, or experience seems to indicate that the very ones hardest 
to avoid are actually the ones having the least noticeable effect, 
and vice versa. 

Under normal operating conditions it was our practice to discard 
a foil or target if the blackeninz became appreciable, This was 
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especially necessary for the foil; evidently the greater diffuseness 
of the beam on the targets and their preater ability to conduct the 
heat away prevented rapid "carbonization" of the targets. 

In connection with the foils, this carbonizing added small 
amounts to the stopping. The effect was checked by determining 
energy loss throush the foils at intervals throughout their use, 
Figure 10, which gives the energy loss in the various foils as 
function of the entefhz proton energy, shows this effect clearly. 
Note especially the progressive increase in stopping power by the 
.05 ril copper foil for subsequent days (on the order of 2 kev per 
operating day for 1 Mev protons). Also note the radical shift in 
energy loss through the „05 mil nickel foil (#2), which occurred 
after a rather severe carbonizing due to release of gaseous hydro- 
carbons into the gas chamber from incompletely dried applications 
of glyptal. 

Foil stopping power data were discarded when it was evident that 
such carbonizing had contributed to abnormally high values of enerzy 
loss. This data is needed, however, when the carbonized foils are 
used as windows for the experiments on gas stopping power. It should 
be mentioned that the greatest problem did occmr in the gas stopping 
power experiments, In the foil stopping power experiments, the whole 
system was continuelly pumped by the vacuum pumps; but in the gas 
experiments, the gas chamber was closed off and any hydrocarbon mole- 
cules trapped therein had many chances to strike the hot foil and be 


decomposed am deposited there. The examples of severe carbonizing 
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nentioned in the preceding paragraph occurred during experiments on 


BASCSe 





Accurate Cormutation of Mean Ener of Protons in liediun, 





It is usually qistomary in giving experimental results for 
stopping power to consider that AE/Ax equals dE/dx for a value 
of E which is the arithmetic average of the initial and exit ener- 
gies of the bean in the stopping mediun (= Nm mor 
exanple, in case the stopping material is a foil). This approxi- 
mation is appropriate when ME is small compared to E. Warshaw 
(19/9) has commented on this matter and proposes a farmula to give 
a correction to AE/ Ax so that the correct value for stopping 
can be obtained at the arithmetic mean value of enercy. His formula 
assumes that the stopping power curve over the region from Eg to 
E is a straight line, He finds that the corrections this com- 
puted are less than his probable error but worth including to insure 
a correct shane for his curves, 

Our values of both AE and E are in most cases larger than 
Warshaw's, and we may also find a slisht correction advantaceous. 
Rather than Warshaw's formula correcting stopping puwer ue prefer one 
which provides a correction to the arithmetic mean enersy to give the 
true nean (designated as Eg) . In deriving the formula for this 
correction let us make an assumption as to the shape of the stopping 
power curve which is more accurate than the linear shane assumed by 
Warshaw. Assune that 


(IV-A-3 ) 
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Theory indicates, as we have noted before, that the exponent Y 
should have a value of 0.5 (See Section III-G). Experimental re- 
sults (Madsen and Venkateswarlu, 1948a; Kahn, 1953) show that Y is 
slightly dependent on Z and on E, but has in the cases wherein 
we are required to use it (E 600 kev) an average value of about 
0.45 . 

The thickness of the medium traversed, Ax , is easily seen 


to be given by: 


1 
25 dd Eo 
T Y 
Ax = - 5 E - dE 3 
dE/dx k (IV-A-4) 
E, Bo 
1 UOI Y + | 


Now the correct value of mean energy, En › to use is that for which 


dE E? - E 
NEU = AEAw « 22 š (IV-A-6) 
Ax 
Solving for Ax and substituting equation (IV-A-3) gives: 
tr s 
Ax ΤῈΣ (IV-A-7) 
k E. 


If we equate the two expressions for Ax and solve for Ei» we obtain: 


кд = y 4- | l/y 
Б = m OPE с o cm (IV-A-8) 
m Y + 1 ES = Eo 


If we replace Ef by E, + SAE, and Eo by F, ` LAE, where 


E, is the arithmetic average, we find: 





8 
Yel mp / 
€ a dam ^ 4. Y 
EQ = l ` (Es + ΔΕ) (Ea SAT) : (IV-A 388 
T AE 
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If we expand the parenthetical expressions sccording to the 


binomial formula, we find: 


l^. 


m , 


2 
E = E, + 1 Y ( Y - 1) (A E)” + other terms 
24 pe-Y A A "a 


š y (IV-A -10) 
where the other terms are successive terms in (AE) /EH-" , 


for n =æ 2, 4, 6, etc., with fractional coefficients for the terms. 
For ME appreciably less than E , this series converges repidly, so 
that the other terms in the above equation may be neglected. 

Then let us expand ene bracketed expression in equation (IV-A-10) 


once again by the binomial formula, to obtain: 


А І 2 
RE: - Ü PES + additional terms , (EV-A-11) 


where the additional terms can once more be dropped when AE is 
appreciably less than E , 


Then we see, substituting 0.45 for Y, that 


E = E - Š E , (IV-A-12) 
where à 
~ -55 (AE) 
dE = E MESE : (IV-A-13) 
24 E, 


We use this correction and apply it below in cases where it 
exceeds 1 kev. 
5. Corrections Required for Impurities in the Materials Studied. _ 
One possible source of consistent error which should be care- 
fully checked is the impurities in the elements studied. 


The metals are quickly disposed of. A semi~quantitative 


-A 





spectrosraphic malysis was carried out with sample foils fron 
the same set as those used, with the following results: 

(a) For copper, the impurities present are as follows: calciun 
with the order of magnitude of 0.1%, magnesium less than 0.1%, and 
a trace of silicon. 

(b) For nickel, the impurities present are as follows: cobalt to 
an amount between 0.1 to 0.5%, plus traces of magnesium, molybdenum, 
beryllium, and iron, 

The error in stopping power is not as great as the percentage of 
impurities. Each of the impurities tends to produce an error only 
to the extent that its stoppin: power varies from the stopping power 
of the prime constituent, so that the closer the impurity is to the 
primary element in the table of elements, the less the effect of the 
impurity on the stopping power measured. With these considerations 
in mind, it becomes clear that the impurities in foils of either 
nickel or copper are such in quantity and type that no appreciable 
error can result by considering tne foils as pure, ; 

The gases hovyever require nore careful study. The gases” as 
they are supplied by the company are guarancvecd to be pure to the 


following extent: 


1 
The analysis was carried out by lr. Sam Wohifort at the Spectro- 
graphic Laboratory, Department of Chemistry, Ohio State Universi- 
ty, through the courtesy of Professor W.li. Laclievin. 

2 


The nitrozen and arzon were obtained in pressurized K cylinders 
from the Linde Air Products Company; while the rarer gases cane 
from the same company in glass flasks at approximately atmo- 
spheric pressure. 
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(a) Commercial nitrogen used is stated to be 99.7% pure. 
(b) The argon supplied is stated to be 99.6% pure. 

(c) The neon impurities are negligible, 

(4) The krypton may have as much as 1% xenon present. 
(e) The xenon may have as much as 2% krypton present. 

Since the above amounts are maximm values, the exact extent 
of the impurities for our samples was still somewhat in doubt. 
There existed a still nore serious factor in connection with the 
latter three gases, which were supplied in one liter or one~half 
liter glass flasks, Difficulties were encountered in the attach- 
ment of the flasks to the system in such a way as to be absolutely 
air-tight when the duct system is under partial or camplete vacuum, 
The results were that a certain amount of air contamination was 
probable for the gases in each of the three flasks, before final 
pas-tLishtness was attained. It was considered necessary therefore 
that an analysis be nade on the gases in the flasks after the pos- 
sible contamination occurred, | 

The analysis was carried out пу nass-spectrograph” » and results 
are as given in Table XVI, 

In order to determine the effects of these impurities, it is 
necessary to note that the amount of matter traversed is indicated 
by the measured values of path length, gas pressure, and gas tempera- 


ture. Since the pressures used are quite low (a few centimeters of 


The analysis was carried out by Ir. Leonard Mak, at the Mass- 
Spectrogranhic Laboratory, Department of Chemistry, Ohio State 
University, through the courtesy of Professor H.L.Johnston and 
Dr. David White, 


mercury) the ideal gas relationships are valid, and the variables 
cited above are sufficient to give the mmber of molecules per square 
centimeter of cross-section perpendicular to the proton path. For 
the rare gases, the atonic stopping power is the same as the molecu= 
lar stooping power; and since the impurities are given as volune 
percent (equivalent to mole percent), the corrective factor to be 
applied to the experinental results is the relative molecular stopping 
power of the impure gas as compared to that for the pure gas. 

Table XVII gives the values assumed far atomic and molecular 
stopping power for the various elements and compounds concerned. It 
has been necessary in some cases to make rough preliminary estimates 
of stopping power fron the data taken in our experiments, so that 
the values taken can be justified only by the final answer.  Fortu- 
nately, hizhly accurate values are not necessary in this particular 
set of computations. The sources of data aro indicated in the table. 
| Table XVIII tabulates the calculations of the corrective factors 
for the gases neon and xenon. It is quite evident fron the analysis 
of the gases in Table XVI that no appreciable leakage into the systen 
occurred for the lrypton experiments, so that no correction is re- 
quired for this gas. The correction factors computed rust divide the 
experimental results to obtain final cprected results. 

Ме have neglected any erfect of possible impurities in nitrogen 
and argon, This is justified for two reasons. First, the gas cham- 
ber vas supplied with nitrogon or argon directly from the pressurized 


tank, so that the filling systen up to the inlet valve of the chamber 
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Table XVI. 
Impurities in Volume Percent in Gases, As 


Determined by Mass Spectrograph 


Primary Gas 
Impurity 

H50 2.72% - 0.02% 
05 Dol .011 Murs 
Α „ОО - ο 
соь «09 2005 - 
Kr 7219 ~ 1758 
Ne € .026 . 07 


Table xVII. 


Approximate Values of Atomic and J:olecular Stoorins 
Power of Various Gaseous Constituents for Protons at 
snergies of 500 and 1000 kev 
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Atom or Atomic Molecular 
yolecule Stopring lower Stopping Power 
( еу-сп& =< 107°) (ev-cn“ x 1019) 






o00 kev | 1000 ke 





ne 159 (5) 

бо - - (2) 

6.0 16.0 9.8 (1) 
6.8 17.6 1150 (5) 
10.1 1000 6.45 (3) 
14.4 ТАА 9.8 (1) 
2n 1 21.1 14.2 (1) 
ee З 04.2 16.25 (1) 
- - 2,0 6.7 (4) 

- В 24.85 15.35 (4) 





on 


(1) Estimated from „reliminary results, present research. 

(2) Extrapolation from data given by Livingston and Bethe (1937). 

(3) Interpolated value, zccording to atomic number, from values 
obtained by references (1) anud (2). 

(4) Cbtained by addition of atomic stop ing powers. 
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Table XVIIL. 
Computation of Impurity Correction Factors, to Convert Experimenta!ly 
Determined .tomic Stopping Fower to the Correct Velue - Velid for 


Trotons of 500 and 1000 kev. 





р + Pa Qa e 


1 e ΤΊ. PME 
l'olecule or (103) Relative 97°, |Vol.| (Vol.%) x (Rel. 


: 
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(a) Neon 
Heon TOS]. 5205 1,00 1,00 86,15 86.15 
HoO 11,0 6.7 1.09 1.04 Doy 2,94 
No 16.0 9.8 1.53 1.02 cT 123,9 
Oo 17.0 11 1,74 1.70 dad 1,91 
IN 14.4 J393 1.43 LE . L „14 
COs 24.85 19235 2.46 2456 “5 1.48 
Kr ol 14.2 2.09 2.20 1.2 2.91 
100.0 108.03 
(b) Xenon 
Xenon 24.2 16,25 1.00 1.00 95.3 95.5 
N 10.1 6.45 „42 „40 21 . 04 
No 19:0 uoo . 06 . 60 >23 2,44 
On 17.6 11,0 ele 208 vy SOL 
Kr 21.1 1402 .87 το; me . 17 
100.0 95.46 
Summary: To correct exverinental data, divide by factors: 
Prime Proton Lnergy 
Stopping 500 kev 1000 kev 
Substance __ 
Neon 1.080 T2079 
Xenon . 982 





was under a positive pressure relative to atmospheric. This pre- 
vented inward leaking of air, so that most of the impurities pre- 
sent could be attributed to those originally present in the gas 
containers purchased. In regard to these original impurities, we 
note that the percentase allowed is miall, Furthermore such in- 
purities are primarily air constituents, Table XVII shows that 

the molecular stopping power of nitrogen, oxygen, and arzon are all 
very close together, so that mall amounts of one gas in the other 
as impurity would have an extremely small effect on stopping pover. 
б. Possible Errors Due to Air Leakare Durinc Experiment, 

It was discovered durins most of the exoerinents on the gases 
that the gas chamber pressure increesed slowly, by varying amounts 
at different times. At the worst, this change amounted to about 254 
of the initial pressure. This indicated a leak of sone kind into 
the gas system during operation, Continual monitoring of the pres- 
sure indicated that such leaking occurred at a fairly constant rate 
for a single filling, although the rate varied from one day to the 
next. The errors due to such leaking could be allowed for therefore, 
provided the type of sas leaking were mom, It might be possible 
that the leakages were either of air, or of the sas being studied, 
through the fillinz system. The possibility also existed that a 
combination of both effects in unlmowm proportions occurred. 

Great precautions were taken to reduce this effect to a ninimum, 
and frequent observation was made on the rate of leakage into the 
chamber with no gas present. Except for a couple of times when 


definite air-leakgzo vas showm (data taken on these days showed radi- 
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cal departures from other data on the same gas, and were discarded), 
no appreciable pressure change was discovered in the evacuated and 
closed gas chanber over periods of several hours duration. 

Despite these checks for air leaks, tho E fect persisted during 
periods that gases bèing studied were in the system. After sone dis- 
cussion with other persons in the departmentsof physics and chenistry 
who worked with high vacuuns, it was concluded that the effect was 
due to the entrapment of minte bubbles of sas in the stop=cock 
grease at the stop-cock which occurred during the filling process. 
Such bubbles sradually traveled to the low pressure port and dis- 
charzed the gas into the chamber. The prominence of such effect de- 
pends upon many factors - type of oe used, type o? vacuum 
grease used, amount of vacuum grease applied, pressure with which the 
valve of the stop-cock is forced onto its seat, It is conclvded that 
the nain difficulty was due to use of a type of srease used for ordi- 
nary chemical applications, but not especially siited for high vacuum 
work, For future workers, Apiezon, Type П, grease, apolied strictly 
in accordance with instructions on the label and in sparing quantities, 
is recormended, 

Notwithstanding the failure to reduce the effect of pressure 
change to zero, the experiments made and precautions taken insured 
that a nezlizible proportion vas due to air leakage. All that is 
necessary for our analysis is to know the precise value of gas pres- 
sure in the chamber at the tine that a resonance peak is measured. 
Since the pressure in the system at various times during the se- 


quence of runs was recorded, and followed very closely a linear 
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change with time, pressures and variables dependent upon time can be 
determined through a simple linear formula for each gas filling. 


These formulas will be given in the compilation of data where they 


apply. 


B. Stopping Power Resuats - Foils, 


1, Determination of Inercy of Displaced Peaks, Foils Inserted. 
Tables XIX and XX show the tabulated data and results of deter- 


mination of displaced peak energy values (E, ). The values of I' , 
the displaced peak position on a magnet current scale, are given by 
the yield plots (Figure 8(b), for exemple). The values of E,' are 
bhen deternined fron the calibration curves (Figure 9). The proper 
calibration curve to use on any particuler day is determined by one 
or more calibration points taken in the same sequence of runs. This 
is discussed in detail in Section III.H, above, 

2, Computation of Foil Stopping Power. 

Tables XXI and XXII indicate in tabulated form the computations 
for stopping power of nickel and copper, respectively. The defini- 
tions and methods of computing the variables specified in the colum 
headings are already completely set forth. Stopping power is com 
puted both in the "experimental" units and also, as atomic stopping 
power, in those units closer to basic theory. A column is provided 
for tabulation of values of 0% 2!/. 2. Since Bohr's theory for 
stopping of slov particles in hoavy nedia provides thot stopping 
power is proportional to α1/ 2 (Equation I-B-74), such a tabulation 
gives a better means of comparing the results of stopping power ex- 
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Table XIX. 
Determination of Shifted Resonance Pesk Positions in the Energy 
Spectrum, with Various Nickel Foils Inserted into the Proton Beam. 


ee O A -αᾱ-. -— RNA mas. sma am. Q --- «αν -ᾱν- 


Shifted Peak Position 


> * 
Foil Date Eo I’ log I*| “log ES m 
(Lev) (amp) (Мету) 





„О5п11 Ni-l 8-29-52 „6714 .7348 1.8948 1.9241 .8397 
" ú .8756 .8696 1.9395 0.0115 1.0268 

" 9-26-52 «9875 .9150 1.9614 0.0517 1.1264 
.O5mil Ni-2 1-19-53 .6696 7908 1.8981 1.9245 „8404 
4 " .8740 8776 1.9435 0.0107 1.0249 

E D „9558 «9015 1.9550 0.0330 1.0790 

: " 1-21-53 „6696 «7955 1.9006 1.9274 .8461 
E n " „8740 «8805 1.9456 0.0130 1.0504 
ΚΕ т " „9358 .9075 1.9578 0.0360 1.0864 
.04 mil Ni 10-16-52 „6714 .7720 1.8876 1.9084 «8098 
n п „8756 .8620 1,9255 0.0027 1.0062 

i я „9574 .8860 1.9474 0.0261 1.0619 

n 10-17-52 .6319 .7560 1.8785 1,8905 .7771 

" " . 9875 .9069 1.9576 0.0455 1.1053 

n 10-23-52  .4441 6681 1.8248 1.7870 «6125 

m r 10-28-52 „9374 .8868 1.9478 0.0270 1.0641 
Л т 10-31-52 „9574 .8870 1.9479 0.0283 1.0673 
.04 mil Ni 12-19-52  .9875 .9145 1.9612 0.0440 1.1066 
π 1-9-55 „6696 7765 1.8901 1.9080 .8091 

" " „8740 .8650 1.9370 1.9986 . 9968 

.03 mil Ni  1-22-53  .6696 „7519 .8762 1.8840 ‚7656 


1 
" " „8740 „427 1.9257 1.9805 .9561 
" " ‚9558 „8707 1.9399 0.0079 1.0184 


чире’ ре τα a wm e 





a mm 


Notes: * The dates not only serve to identify the various sets of 
data extracted from vlotted yield curves, but also 
indicate the proper calibration curve to use (Fig.9). 


# Obtained from calibration curves (Fig. 9), with use of 
Tables XIV and XV. 


j- Appreciable carbon accumulation indicated. 
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Table XX. 
Determination of Shifted Resonance Peak Positions in the Fnergy 


Spectrum, with Various Cooper Foils Inserted in the Proton Beam. 


Shifted Peak Position 


Foil Date Eo I! log 1* | log ES 
(Hev) (amp) 





.05 mil Cu  12-2-52 „5507 .6310 1.8000 1.7372 „5460 
" " .6763 «7860 1.8954 1.9215 .8346 
: D «8799 .8735 1.9413 0.0080 1.0186 
Г f .9415 «8970 - 1.9528 0.0298 1.0710 
п 12-35-59 «6919 .7705 1.8868 1.8996 .7956 
f i «9075 .9212 1.9644 0.0485 1,1182 
п 12-18-52 „4419 «676 1.8299 1.7956 .6246 
Bo" " .6696 «78584 1.8940 1.9188 .8?95 
£ " £ .8740 .8710 1.9400 0.0073 1.0170 
#4 i .9558 „8960 1.952535 0.0307 1.0732 
gon 1-12-53 «6696 .7860 1.8954 1.9205 98323 
M t " (pn) .6696  ,7855 1.8951 1.9210 .8337 
O.lrilCu 1-22-53 .4419 .7580 „8797 1.8877 .7721 


8 
9287 1.9819 • 9592 
9 


І 
" " .5696 .8485 L 
1.9644 0.0515 1.1259 


" | n .8735 9514 
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# Appreciable cerbon accumulation indicated. 


See also the notes in Teble XIX. 
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# Appreciable carbon accumulation indicated, 
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Table YXI (Two Pages) 
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0.640 


tt 
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0.977 
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671,4 
875.6 
987.5 


874.0 
935.8 
669.6 
874.0 
935.8 


671.4 
875.5 
957.4 
631.9 
987.5 
444.1 
957,4 
957,4 


987.5 
669.6 
874.0 


669,6 
874.0 
935.8 





1086,4 


809.8 
1006.2 
1061.9 

777 ot 
1105.3 

612.3 
1064.1 
1067.5 


1106.6 
809.1 
996.8 


765.6 
956.1 
1018,4 
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Stopping Power Computation for Nickel 


156.4 
150.6 


138.4 
130.6 
124.5 
145.2 
117,3 
168.2 
126.7 
129,9 


119.1 
139.5 
12248 


96.0 
82.1 
82.6 
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— -ρ. 
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Ea 


(kev) 


795.6 
951.2 
1057,0 


755.0 
949,4 
1007.4 
757,8 
906. 2 
1011,1 


740. 6 
940.9 
90886 
704,5 
1046.4 
528.2 
1000.8 
1002.4 


1047.0 
739,4 
935.4 


717.6 
915,0 
977.1 











Table XXI. (Concl.) 


Stopping Power Computation for Nickel 


Б. -Ὢ << u = - — eo — -- 














Stovping O” , АФ. Prob. 
Power St pg . Power 97 21/5 Error, 
(еу-ст2) σ 


(kev-en*/mg) 
















A A Ó — = — e 


757 146.3 14.27 4.70 1 
951 131.5 12.51 4.12 2 
1057 120.7 11.77 5.88 2 
755 149,3 14.18 4.67 1 
949 128.3 12,91 4.12 2 
1007 LA 11,88 5,91 2 


Carbon accumulation renders invelid. 


$1 " " τ 


n τ tt MU 


10 negl. 741 141.7 15.80 4.55 249 
11 " 941 15527 13.02 4.29 ο 
12 " 1000 127.4 12,42 4.09 2.0 
13 E 704 148,7 14,48 4.77 "eO 
14 i 1046 120,4 11.74 Ode QD 
15 1,2 527 172.0 10% 8 or o2 wS 
16 Carbon accumulation renders invalid. 

17 τῇ T 1 1 


ngl. 1047 122,0 11,89 3.92 2.5 
759 142.8 1,91 4.58 2.0 
935 125.7 Tosco 4,04 2.0 
718 150,0 14.62 4.81 520 
915 128,2 12.50 4,12 3.0 
Ө 129.0 12.57 4,14 3.0 


— e v <“ o < = €———— 
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4 Values tabulated in these columns to be multinlied by 10 


7, 


—f Y — 


Line 
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Table XXII. (Two Pages) 


Stopping Power Computation for Copper 


Foil 
Number 


209 nir cu 
71 


0.1 míl Cu 


e" 


tt 






Foil 
Thickness 


(mg/em”) 


2.20 
" 


350.7 
676.3 
879.9 
941.5 
631.9 
987.5 
441.9 
669.6 
874.0 
935.8 
669.6 
669.6 


441.9 
. 669.6 
873,5 





Eo ES 
(kev) (kev) 


546.0 
834.6 
1018.6 
1071.0 
793.6 
1118.2 
624.6 
829.5 
1017.0 
1073.2 
832.3 
833.7 


772.1 
959.2 
1125.9 


# Appreciable carbon accumulation indicated. 


-/ος 9 


AE 
(kev) 


195.3 
158.3 
138.7 
129.5 
161.7 
130.7 
182.7 
159.9 
143.0 
137,4 
162.7 
164.1 


350.2 
289.6 
252.4 


(kev) 


448,4 
795.5 
949.2 
1006.2 
712.8 
1049.6 
533.2 
749.6 
945,5 
1004.5 
751.0 
751.6 


607.0 
814.4 
997.8 


Table XXII. (Concl.) 


Stopping Power Computation for Copper 


ae ey aa «αν Ss 

























Stopping # 0” , At. # 
Power St pg . Power o7/ 71/5 
(kev-em^/ng)|  ( ev-cm^ ) 
1 169.0 17.86 5.815 
2 negl. 755 157.0 14,48 4.715 
à " 949 120,0 12.68 4.13 
4 б 1006 112.0 11.83 5585 
D n 713 139.9 14.78 4.81 
6 H 1050 713.1 11295 3.89 
7 1.4 532 19820 16.70 5.44 
8 Carbon accumulation renders invalid, 
9 " n 
10 n " 
11 tt tt 
12 η τ! 
15001 15.87 5.165 
131,6 13.90 4.525 
114.7 111 3.94 








D areas: > ee eee cc cc c ccc c cci eee i edi d 


# Values tabulated in these columns to be multivlied by 10719, 
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periments on different elements. 
3, Discussion of Errors in Experimental Data and Results, 

Probable errors in the individual values of stopping power ob- 
tained are difficult to judge Secum on. but even an approxinate 
idea of such deviations would be useful. | 

Values obtained fron the literature for the resonance energy 
values used are probably correct to about # 1 kev on the average. 
Some are known with better accuracy, since they are quoted to tenths 
of a kev, Many other peak values are determined in relation to 
those few peaks whose energy have been accurately obtained by abso- 
lute methods, and thus their error may be mich larger. ‘the value of 
l kev is considered a rough average for all those used herein, 

Errors in target thickness are likely to be of the same order 
of magnitude. However, although such errors affect the values taken 
for Eg , they cause at the same time an erroneous shift in the cali- 
bration curve used which introduces an error of the same order of 
magnitude in the value of Eg', Thus for AE , these two errors 
largely cancel, We will therefore neglect such errors. 

The important errors to be considered in E,' are due to error 
in obtaining the proper value of I' (value of shifted peak position 
in the magnet current soectrum) and to an error in choice of the cali- 
“bration curve used. In regard to the former it must be noted that 
the stopping phenomenon not only shifts the peak but also broadens 
it, so that the position of the maximum point is more difficult to 
estimate. Whereas the undeviated peak can be judged to an accuracy 


of about 0.1 milliameres, or about 0.2 kev, the deviated peak has a 


Ef ) 
AR > 





greater probable error = probably on the order of 1 kev. 

| In regard to the error from calibration inaccuracies, let us 
refer to Section JII.H wherein a typical value of probable error 
is about 1.7 kev by use of the mean curve. We can improve on this, 
however, since we determine one or more calibration points during a 
given day's sequence of operations. A value of 1 kev should be 
about the proper value for this error. 

Summarizing bat we seo that the probable error in EJ! 
is about 1.4 kev anc about 1 kev for Eg, so that the probable 
error from the difference of these two is about 1.7 kev. We fcel 
it necessary to add to this an estimate for surface impurities, 
primabily carbon, of about 1.4 kev (see Section IV.A.3.). This 
gives a final probable error in ME of about 2.2 kev. 

The percentage error which this figure provides depends on the 
value of AE, Reference to Tables XXI and XXII indicates that the 
probable error varies from 0.7% for thick foils and low proton ener- 
gies to 2.7%-for thin foils and hich proton energies. to take into 
account the possible increase in path length at low energies, let 
us raise the lower figure to 1 lev, 

It has already been shown in Section III.B. that the probable 
error in foil thicknesses is about 1.1%. This means that the 
probable error in values for fonl stopping power veries between 1.5% 
and 3%. Tables XXI and XXII tabulate the estimates of probable 
error for individual values of stopping pover. 

It must be emphasized that the probable errors given are based 


on the assumption of all errors conforming to the Gaussian theory of 


4 
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errors, with positive and negative errors equally likely. Further- 
more this analysis does not guarantee that unknown errors of either 
a consistent or random nature have not caused the tabulated values 


of probable error to be less than they should be, 


C. _Stoppinz Power Results - Gases. 





1. Έπεγσ Loss in Window Foils. 

It is evident that in obtaining the energy loss throuch the gas 
in the chamber, it is necessary not only to subtract E, from T 5 
but from this difference to subtract the energy loss in the window 
foil (Equation II-B-9). Since for windows we used the same foils 
on which the stopping power measurements were made, we can determine 
the value of window foil energy loss, Alp , by adis. use of a 
chart whereon is drawn for each foil used a curve showing foil 
energy loss versus proton emergy at the foil entrance, E,'. Using 
the data provided in Table XXI and XXII, this can be readily accom- 
plished. Figure 10 gives these curves. It should be noted that 
certain points are plotted which indicate clearly the effect of 
carbon accumulation. Such data enable one to make a better estimate 
of А Ер so as to include carbonizing when appropriate. 
2, Energy of Displaced Peaks, Gases and Window Foils Inserted. 

These values are determined precisely as they are for the case 
of foils alone (Section 1V.B.1). Tables XXIII-XXVII, inclusive, 
show this determination, 


3, Determination of the Formmla for Atomic Stopping Power of Gases. 





The fundamental conversion formulas for stopping power have al- 


—/ wI 





rcady been evaluated. It is necessary however to develop a formula 
which will give atomic stopping power as a function of gas tenpera- 
tore (T), pressure (p), and beam energy loss. ‘The other factor, 
path distance, is already known from measurement of the length of 
the gas chamber (Fig.ıre 4). This diStance, A X, is 10.56 cm. From 


Section I.A., we see that 


3 ᾿ 
Q- (ev-en* = 10 . BE (IV-C-1) 
( ) — 10.56 3 


where AE is measured in kev, 
The ras pressures used are so low that the ideal gas relations 
are sufficiently accurate for our purposes. Then, under standard 


conditions, 


N = {ο x 1042 . n (τη 6”) 
P EN j 


where n is the number of atoms per molecule of the gas considered, 


Furthermore, 


| е 2 s 
N = dE e ` ose i (1V-c-3) 


where p is measured in centimeters of mercury, and T is measured 
in degrees Kelvin. 


Combination of the above formula gives the relation sought: 
c = СЛЕ 9 


where (IV-C-4) 
G = T 
1,021 C l "Dan 


Since p and T are functions of time (Section IV.A.C), it is possible 


to determine G as a function of time for each gas filling. Table 
A [q - 





Table XXIII. 


Determination of Shifted Resonance Peak Positions in the 
Energy Spectrum, with Nitrogen Gas and !!indow Foil Inserted 


in the Bean. 


— e r - m dili -. ee a s ee 


Shifted Peak Position 


E 
Foil Date (Mev) It log 1'] log ES ES 
(amp) (Mev) 


Window 





.04 mil Ni 10-57-55 «5458 6820 338 1.8040 „6568 


1.8 
" " ¿4441  .7237 1.8596 1.8545 «7150 
P 7 «6006 .7846 1.8946 1.9225 «8566 
p | .6714 .81з2 1.9102 1.9525 .8964 
" " «8756 .8936 1.9511 0.0318 1.0760 
D S .9374 „9164 1.9621 0.0535 1.1306 
.04 mil Ni 10-28-52 .3438 .6852 1.8358 1.8085 „6454 
| Л „4441 7255 1.8606 1,8565 „7183 
" " .6006 „7856 1.8952 1.9234 «8585 
" " 6714 „8150 1.9112 1.9544 «9005 
" E .8756 «8940 1.9513 0.0321 1.0767 
à | .9574 — .9168 1.9625 0.0535 1.1511 


.08 mil Ni 1-23-53 9413 .7050 1.8482 1,8295 .6753 
" " .4419 .7413 1.8700 1.8716 . 7440 
" " «6696 | .8269 1.9175 1.9622 „9166 
" " .8740 «9045 1.9564 0.0368 1.0884 
" " „9358 „9279 1.9575 0.0580 1.1429 





TSE SS i ee EE See SE eee Ee ν΄ ШР СӘ RR o 


Obtained from calibration curves (Fig. 9), with the use of 
Tables XIV and XV. 
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Table XXIV. 


Determination of Shifted Resonance Peak Positions in the 
Energy Spectrum, with Argon Gas and Window Foil Inserted 


in the Beam. 





чч - - LU Ñ 





т-а» πα р ар ee (> δι αμ. "πι v ee 


Shifted Peak Position 








Wind ow 
Foil 












.04 mil Ni 10-29-52 
" " .4441 


„04 mil Ni 10-30-52 „5458 
tt tt 4441 
1 t" š 6006 
" " „8756 
п п „9574 „9280 675 0.0658 


pa! 
© 
C 
CO 
j=- 


.04 mil Ni 10-30-52 .8756 . 9080 0.0474 - 


p.m, 


e 






του ο κο.  - a 9 ъч u 





ασ ο» - ж зш A O A AA TA ы эв 


* obtained from calibration curves (Fig. 9), with use of Tables 
XIV and Ve 
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Determination of Shifted Resonance Peak Positions in the 


Table XXV. 


< 


Energy Spectrum, with Neon Gas and Window Foll Inserted in 


the Beam. 


Shifted Peak Position 


Window 
Foil Date з, I log I' y BE ES 
(amp) (Mev) 
# ,OSmil Ni-2 1-20-53 5415 .7164 1.8552 1.8453 «7005 
п. " 4419 .7572 1.8792 1.8910 .7780 
à „5986 .8180 1.9128 1.9553 .9022 
" " .6696 „8440 1.9263 1.9810 . 9572 
" " 8740 „9210 1.9643 0.0555 1.1511 
a " . 9358 69452 1.9755 0.0746 1.1874 
# ,O5mil Ni-2 1-21-53 .3413 .7220 1.8585 1.8475 „7059 
п " «4419 «7617 1.8818 1.8921 . 7800 
" " e 6696 „8469 1. 92 78 1 9802 е 9554 
п on «8740 .9250 1.9661 0.0535 1.1311 
" " .9358 «9466 1.9762 0.0725 1.1817 


+ 
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* 
Obtained from calibration curves (Fig. 9), with the use of 
Tables XIV ana XV. 


fr Foil was blackened by heavy carbon deposit during runs on 


20 January. Target thickness was checked before use on 
21 January (Figure 10). 
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Determination of Shifted Resonance Peak Positions in the 


Table XXVI. 


Energy Spectrum, with Krypton Gas and Window Foil Inserted 


Window 
Foil 


" 
" 
1? 


.05 mil Cu 


et 


Date 


1-12-55 
" 





in the Beam. 


E 


«4419 . 7518 


„5986 „8130 
„6696 „8422 
.8740 «3506 
5413 . 7070 
„4419 «7492 
„5986 „8110 
„6696 „8404 
„8740 „9170 
. 9358 «3594 


SSeS Sf P 





кҥн! | 
& 

о «D «o 0 
е 
O 
= 





Shifted Peak Position 


O It 1ος 1" log Ej Es - 
(Mev) Ë 
(amp) (Mev) 


„7617? 


„6811 
«7610 
.8867 
„9495 
1.1228 
1.1702 


И Obtained from calibration curves (Fig. 9), with use of Tables 


XIV and XV. 
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Table XXVII. 


Determination of Shifted Resonance Peak Positions in the 
Energy Spectrun, with Xenon Gas and Window Foil Inserted 


in the Bean. 





С чаннар заро с ары l Б ЫА - 





Window Shifted Peak Position 
Foil Date ἘΝ r | 108 1» | *208 Es E? 
(Mev) (amp) (Lev) 

$05 mil Cu 1-14-53 „3415 .7270 1.8615 1.8572 . 7198 
" " ¿4419 «77286 1.8801 1.9080 „8091 
" " 5986 .8268 1.9174 1.9645 „9215 
" " 26696 „8561 1.9325 1.9932 „9845 
" " ‚8740 „9292 1.9681 0.0613 1.1516 

005 mil Ni-2 1-19-53 5415 .7351 1.8663 1.8635 . 7303 
" " 4109 .7747 1.8891 1.9073 .8078 
" " .5986 .8300 1.9191 1.9647 „9219 
" " .6696 .8581 1.9335 1.9923 . 9824 
" " 8740 .9517 1.9785 0.0791 1.1998 


* 
Obtained from calibration curves (Fig. 9), with use of Tables 
XIV and XV, 
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Table XXVIII. (Two Pages) 


Computation of G(t) , a Factor Required ín Com- 


puting Stopping Power of Gases (Equation IV-C-4). 





* 
Values tabulated are to be multiplíed by 10 


~/32-- 


Tabs. , Hour, t, $ 
Date Initial: Initial: G(t) 
Final : : : | Final 
(a) Nitrogen 
10-27-52 298.5 2.02 7239 - 72.5 
299.5 2005 7829 
10-28-52 297,4 2.076 70.1 9:45 70,1 - 0,847(t-9.25) 
298,9 2.20 66.5 14:00 
1-23-53 298,1 3.915 37,25 - 37.25 
b) Argon 
10-29-52 297.5 529 5961 10:15 55.1 -0.392(t-10.25) 
298.8 5.60 5505 14:50 
10-50-52 29758 3.10 94.0 12:30 94,0 + 0,20(t-12.5) 
A elle 299.2 2909 94.8 16:30 
10-30-52 299,2 ο. Ὁ 90.1 - 90.1 
P.M. 
(ο) Neon 
1-20-53 298.0 4.685 13:27 62.25-0.925(t-13.45) 
298.2 4.815 15:14 
1-21-53 298.0 4,28 11:05 68.2 
298.0 4.39 1125 66.4 
298.0 42515 13:10 64.6 
298.5 4.85 14:00 60.25 
298.5 4.86 14:10 60.1 
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Table XXVIII. (Concl,) 


* 








TOAbs., | p (em), an Hour,t, Й 
Date Initial:| Initial: | Initial: | Initial: G(t) 
Final :| Final : | Final Final 
d ton 
298,5 2.085 109,3 16:45 
1-12-53 298.0 2.097 "142.0 13:00 142.0 - 4.16(t-13.0) 
298,5 2.27 131.6 15:30 
(e) Xenon 
1-14-53 298,5 2.183 133.9 10:30 133.9 - 12(t-10.5) 
299.0 2.482 117.9 11:50 
1-19-55 298,5 2.204 132.7 13:33 132.7-12.84(t-13.55) 
1-19-55 298,4 3.36 87.0 16:00 87.0 


* 


Values tabulated are to be multiplied by 10 


t Ó Ñ: 
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XXVIII shows the values of G(t) for each sequence of runs in 
studying gas stopping power. 
4. Computation of Gas Stonping Power. 

Tables XXIX throu:sh XXXITI, inclusive, give the computations in 
tabulated form for stovpins power of the gases nitrozen, argon, 
neon, krypton, and xenon, The various factors involved have al~ 
ready been fully defined and explained. The results are expressed 
in the same units as the foil stopping power, and with the same end 
in view. (See Section IV.B.2) 

5. Discussion of Errors in Gas Stopping Power Results. 

As in the case with the foíls, we should expect the probable 
error in position of the undeviated peak to be on the order of 1 kev. 
Er-ors in the value of the shifted peak are, as before, about 1.4 kev. 
We have to subtract from tne energy shift the energy loss in the 
window foil. The values of the curves plotted in Fim re 10 are 
correct y 2 kev, and an additional 1.5 kev variation in the 
stopping power of the carbon film is estimated. These probable 
errors combine statistically to sive an over-all probable error in 
AE of about 3 kev. This amounts to about 1$ for large values of 
AE and up to 4% for the smallest values of energy loss in the gas. 

The probable error in pressure readings is estimated to be 1%, 
and the uncertainty in the quality and quantity of inpurities, 
though doubtful, is taken to cause probable errors in stopping 
power of about 1% for nitrogen and argon, and about 2% for the 


rarer cases. 


—4 44 — 


Line Foíl Date 
1 .04 mil Ni 
2 n tt 
G tt 1? 
4 п f 
5 п tt 
& 4: m 
7 .04 mi1 Ni 10-28-52 
8 tt n 
9 п п 
10 "n τ! 
11 p 9 
12 11 п 
13 .03 mil Ni 1-23-53 
14 n 1? 
15 n tt 
1 6 п m 
17 tt 1? 


* 


Table XXIX. (Three Pages) 


Computation of Stopping Power of Nitrogen (3N,) 


Time 


10-27-52 11:00 


12:30 
13:45 
15:00 
16:00 
17:00 


10:30 
11:15 
12:15 
13:00 
13:30 
14:00 


1 
1 


> > P OI 
8392 


1 
1 
1 


ES 


(kev) 





636,8 
715.0 
836.6 
896.4 
1076.0 
1130.6 


643.4 
718.3 
838.3 
900.3 
1076. 7 
1131.1 


675.3 
744,0 
916.6 
1088.4 
1142.9 


ar, 


(kev) 


164.5 
154.5 
141.0 
136.5 
124.0 
120,0 


165.5 
155.5 
142.5 
137.5 
125.0 
121.5 


102.5 
97,5 
86.5 
79.0 
77.0 


E3 


(kev) 





472.3 
560.5 
695.6 
759.9 
952.0 
1010.6 


477.9 
562.8 
695.8 
762.8 
951.7 
1009.6 


573.8 
646.5 
830.1 
1009,4 
1065.9 


Allowance made for slight carbonization of foil on 27 and 
28 October (Figure 10). 


k 
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Table XXIX. (Cont.) 









Correctíon 
Factor for 
Impurities 
(Divide) 





Table XXIX. (Concl.) 


Corrected is OF Atel # o Stopping 
AE Stopping / 72/3 Power 
(kev) | (kev-cm2 /mg) 


оо ο ο 
(л Q@ QI @ O O 


к” с› е OO Ο) - 


8 e 
e ο 
m mm Q ooo 


Oo QO > 
NINAN 


+ 
3 
de 
A 
4 
4. 
2 
5 
a 
+ 
+ 
4 


P. Q Q Q t 
Οσοι 





* Data in this column to be multiplied by 10718 . 


if Data in these columns to be multiplied by 10715 š 
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Table XXX. 


(Two Pages 


Computation of Stopping Power of Argon 


Me (kev) (kev) 

1 .04 mil Ni 10-29-52 10:30 765.1 152.0 613.1 

2 E ч 11:15 839.1 144.5 694.6 

3 " " 12:30 963.8 134.5 829.3 

4 d е 13:00 1017.9 130,5 887.4 

o п п 14:00 1185.8 120.0 1065.8 

6 i 10-30-52 13:15 666.0 164,5 501.5 

7 Ш п 14:00 749.9 154.5 595,4 

8 " " 14:45 868.0 144.0 724.0 

S п п 15:30 1109.2 126.0 983.2 
10 " E 16:00 1163.6 122,5 1041.1 
11 ч " (pn) 17:30 1115.5 126,0 989.3 


Е ИШКА с... 
Allowance is msde for moderate carbonization of the window foil 
(Figure 10). 


Correction 
Faetor for 
Impurities 
(Divide) 


— - з 


1.0 


1 
2 
3 
4 
S 


O O O í m 








Table ХХХ. (Сопс1.) 








i gunt. it 7 1/3 Stopping 
G Stopping Z, Power 
(Multip.)| Power (kev-cm* /mg) 

1 269.3 55,0 14,81 5.66 224 2,0 

2 250.5 54.7 15271 5.24 207 DD 
3 228.7 54.2 12.39 4.73 187 Zoo 
4 216.0 54.0 11.66 - 4.45 176 oO 

5 190,2 53.6 10.20 3.89 154 3.0. 

6 197.7 94.2 14.87 5.68 DIO 

7 151,3 94.3 14.28 5.45 3.5 

8 123.4 94.5 11.67 4,45 4.0 

9 107.6 94.6 10.19 3.89 4.0 
10 103.7 94.7 9.82 3.75 4.0 
11 113.7 90.1 10,24 3.91 4.0 


* Data in this column to be multiplied by 10^ 


f Data in these columns to be multiplied by 10^ 





Pa) 


a 


18 


15 









Table XXXI. (Two Pages) 


Computation of Stopping Pover of Neon 


* 

Line Foll Date Time ES A Ep Ex 
(kev) | (kev) (kev) 

1 .05 mil Ni-2 1-20-53 15:35 700.3 190,0 510.3 
2 и i 13:55 778,0 179.0 599,0 
3 " " 14:15 902.2 165.5 736.7 
4 " " 14:30 957.2 160.0 797.2 
5 п " 14 «50 1151.1 144.5 986,6 
6 " η 15:05 1187.4 140.0 1047.4 
"7 τ 1-21-53 11:05 703.9- 198.0 505,9 
8 " " 11:25 780.0 185,5 594.5 
9 п " 15:10 955,4 164,5 790,9 
10 " " 14:00 1131.1 146.0 985.1 
11 " " 14:10 1181.7 141.0 1040,7 


Р The foil was blackened by heavy carbon deposit during run on 
20 January. Carbonized target checked for stopping power and 
used again 21 January. This effect has been taken into account. 
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Correction 

Line Ea = SE En Factor for 

B(E, + Ez) Impurities 

(kev) (kev) (kev) (Divide) 

1 341.3 169.0 425.8 1,5 424 1.081 
2 441.9 157,1 520.4 Te: 519 1.080 
3 598.6 138.1 667.6 negl. 668 1.078 
4 669.6 127.6 733.4 " 733 1.078 
5 874.0 112.6 930.3 " 930 1.076 
6 935.8 111.6 991.6 " 992 1.075 
7 341.3 164.6 423.6 145 422 1.081 
8 441.9 152.6 518.2 1.0 517 1.080 
9 669.6 121.3 750.2 negl. 730 1.078 
10 874,0 111,1 929.5 " 929 1.076 
11 935,8 104.9 988,2 " 988 1.075 





Table XXXI. (Concl.) 


Stopping 
power 


Stopping 


Power (kev-cr? /mg) 





l 156.2 62.1 9.70 4.50 290 3.5 
2 145.4 61.8 8.99 4.17 269 3.5 
3 128.1 61.5 7.88 3.66 23569 4.0 
+ 118.4 61.3 7.26 3.57 217 4.0 
5 104.7 61.0 6.38 2.97 191 4.0 
6 103.8 60.7 6.30 2.93 188 4.0 
2 162,1 6802 10.37 4.82 310 3.5 
8 141.3 66.4 9.38 4.35 280 3.5 
9 112,6 64.6 7.27 3.38 217 4.0 
10 103.3 60.25 6.22 2.89 186 40 
11 97.6 60.1 5.86 2.75 175 4.0 


* 
Data in this column to be multiplied by 10718 . 


# Deta in these columns to be multiplied by 10-15 , 
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Table XXXII. 


Two Pages) 


Computation of Stopping Power of Krypton 


— s - 





Line Foil Date Time ES AJ е 
(kev) | (kev) (kev) 
1 .04 mil Ni 1-9-53 pn 15:10 761.7 146,0 615.7 
à " " 15:45 885.5 134.5 751.0 
3 в Ы 16:05 947.3 130.0 817.3 
4 " " 16:35 1125.9 118.9 1007.4 
5 405 mil Cu 1-12-55m 13:25 681.1 178.5 502.6 
6 " " 13:55 761.0 169.0 592.0 
7 п " 14:15 886.7 156.0 "730.9 
8 " " 14:45 949,5 152.0 797.5 
9 " " 15:00 1122,8 138.0 984,8 
10 m " 15:20 1176.2 155,0 1043.2 
[ —— | Correction 
Line Hee " Factor for 
Impurities 
(kev) (Divide) 
1 527 1.0 
2 598.6 152.4 674.8 negl. 675 " 
3 669,6 147.7 745,4 M 745 " 
4 874,0 133.4 940.7 " 941 " 
5 341.3 16123 421.9 1.4 420 " 
6 441.9 150.1 516.9 1.0 516 " 
7 598,6 132.1 664.6 negl. 665 " 
8 669.6 127.9 255.5 " 733 " 
9 874.0 110.8 929,4 " 929 " 
10 935.8 107.4 989.5 " 989 " 
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Table XXXII. (Conel.) 






Corrected | * Fo с Ας, Stopring prob. 
Line G Stopping 9721/5 Power Error, 
(liult.)| Power ( kev-cm^ /mg) 4, 


боочо 


Data in this column to be multiplied by 10718 


>* 


Ў Data in these columns to be multiplied by 10719 ° 





Table XXXIII. (Two Pages) 


Computation of Stopping Power of Xenon 








I — 9 араа 


* 





Line Foíl Date Time FS AE, 

(kev) (kev) (kev) 
1 .05 mil Cu 1-14-53 10:35 719.8 183,5 σος 5 
e п " 11:05 809.1 199.5 636.6 
3 ч pl 11:10 921.5 161.0 760.5 
4 " MH 11:20 984,5 154.5 830,0 
5 Ы " 11:35 1151.6 139.5 1012.1 
6 .05mil Ni-2 1-19-53 14:00 730.3 184.5 545.8 
7 ч N 14:15 807,8 174,0 653,8 
8 v М 14:30 OE 9 161.5 760,4 
9 a ч 14:40 982.4 129.9 826.9 
0 i " 16:00 1199.8 156,0 1063.8 








* 








Correctio 





Line ό E Em Factor for 
Impuritíes 
(kev) (kev) (Divide) 

1 341.3 195,0 438.8 2.0 437 
2 441,9 194.7 559.2 1.6 538 
3 598.6 161,9 679.5 negl. 679 
4. 669.6 160.4 749.8 " 750 
5 874,0 158.1 943.0 " 943 
6 341.3 204.5 443.5 ооп 441 
7 441,9 191.9 537.8 1.6 536 
8 598.6 161.8 679.5 negl. 679 
9 669.6 πο Ὁ 748.2 " 748 
10 874.0 189.8 968,9 " 969 
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Table XXXIII. (Concl.) 


H 
if C, At. Stopping Prob. 


Stopping Power Error, 
Power (kev-em^ /mg) 4, 


Corrected 
Line 





1 132.9 26.3 6.96 120.5 5.5 
2 126.9 2931 6.64 115 Seo 
3 125.9 20672 5.48 95 575 
4 123.9 20,22 5.35 92.7 3.5 
5 120.9 16.99 4,49 77.8 4.0 
6 208 126.9 26.4 6.98 121 3.0 
7 194.9 123.7 24.1 6.38 110.5 Sao 
8 164.5 120.5 19.85 5.55 91 З.Ә 
9 160.1 118.3 18.95 5.01 86.9 3.5 
10 193.3 87.0 16.8 4,44 77 =. 
М Data in this column to be multiplied by 10728 9 
# Data in these columns to be multi»lied by 10719 s 
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Thus, the probable error amounts to 2$ to 4.5% for nitrogen and 
argon, and 2,5% to 4.5% for the other sascs studied. The hicher 
values are valid for smaller values of energy loss, and vice versa. 
Tables XXIX through XXXTIT include colums for the probable errors 


computed in this fashion. 


D, Stragglinz in Stopping Power. 
1, Computation. 

We shall consider (%/ Ax as the neasure of stopping power 
straggling, and assume it to be valid for E, , the mean value of 
energy for the protons passing throush the stopping medium, Results 
were considered only for the roils. Results for the gases are not 
easily obtained by using the technique employed for ovr stopping 
power experiments. There are two reasons for this: 

(1) Stragsling results for window foils must be known with great 
accuracy in order to subtract such effects from the overall results 
provided for gas plus window foil. Such accuracy for the foils has 
not been obtained, 

(2) The statistical spread of the resonance peaks, with both gas 
and window foil inserted, is so large that adjacent peaks interfere 
somewhat with one another. This introduces uncertainty in how mich 
of the yield data are ascribable to each individual peak, 

Table XXXIV gives the data and computations for straggling effect 
in proton stopping power. The colum headings have already been de- 
fined and discussed in nrevious sections. The basic formula for de- 


ternination of SE is equation (II-B-10). 


ж / o 





For reasons which have been explained in Section III.G., the 
values of Ns; are not usually large enough to consider. Thus 
the values of nS ο 1,2 are not considered in the final 
determination of 0 , 

In obtaining the values of nn and f from the yield 
curves the following rules of procedure were used (sce Fisures 8 
(a) and (b) for illustration): 

(a) The non-resonance portion of-the yield was estimated and sub- 
tracted from the yield values to give only the resonance portion 
of the yield for analysis. This is estimated as the level which 
the yicld curve seems to approach on either side. Since displace- 
ment of the noneresonance portion of the curve should not affect 
its height, the same non-resonance datum level should be taken for 
corresponding shifted and unshifted peaks. In case the non-reso- 
nance base lines do not appear equal, some compromise has to be 
reached. 

(b) If the peaks appear Gaussian in shape, ve can estimate the 
standard dev ation ovite accurately by taking a half of the peak 
width at 60.6% of its height above the datum level. This is easily 
proved by reference to any table of the error function, or nornal 
distribution. 

(ο) If the peak is not Gaus-ian in shapc, the standard deviation 
must be computed by numerical analysis in accordance with basic 
definition, 


(d) The in®luences of adjacent peake must be subtracted if they 
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Date 


8-29-52 
10-16-52 


τ" 


10-17-52 
tr 


10-23-52 
l- 9-55 


12- 2-52 


t? 


839.7 
809,8 
1006.2 
777.1 
1105.3 
612.3 
809.1 


546.0 
834.6 
1018.6 
1071.0 
829.5 
1017.0 
1073.2 


Computation of Straggling in Stopping 


Table XXXIV. (Two Pages) 


Power of Copper and Nickel Foils 


Foil 





eOSmil Ni-1 1.150 


004 mil Ni 


.05 mil Cu 


n 


0, 977 
tt 


1.155 


Ac 
1.18 


τι 
п 


1.155 


Resonance 
Peak Ener. 


(kev) 


669 


441.1 
669 


340.4 
669 

873.5 
935.3 
669 

873.5 
935.3 


MAA AT AM memo iue “e — 


τ! 
(amp) 


«7848 
«7760 
„8620 
. 7560 
. 9069 
«6681. 
„7765 


„6510 
ο 2960 
«8735 
.8970 
. 7834 
.3710 
.8960 





to thickness because of carbonization. 


| 








JOUAN 
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12 
13 
14 





62.9 
60.8 
92.75 
76.6 
114.3 
87,4 
66.7 


40.2 
66.3 


113.8 


75.0 
77.6 
88.5 


Table XXXIV. (Concl.) 


1,61 
1.56 
1,93 
1.49 
2.13 
1.18 
1.55 


1.05 


1.96 
2.06 
1.59 
1.96 
2.06 


2.59 
2.43 
3.72 
2.22 
4.54 
1.59 
2.40 
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seem to be present to any extent. 

For the sake of comparison with theory, values of straggling are 
plotted in Figure 15. 
me Discussion of imrors. 

It is difficult to make good estimates of probable errors of 
straggling results, since many factors of indeterminable size are 
present. The best estimate of the size of the random errors can be 
gained by noting the spread in the results as tabulated and as 
" plotted in Figure 15. It will be noted that the values for copper 
are rmch more self-consistent that those for nickel. This may be 
due to the fact that the data for copper were in seneral talen at a 
later period than the data for nickel, and the improvement in tech- 
nique and procedures obtained with >ractice tended to make later 
results more precise. Оп the other hand, the amount of data is not 
sufficient to insure that such self-consistency for the copper re- 
sults is more than pure chance. 

Factors leading to larze errors in the strarglins results include 
such matters as are listed: 

(a) There is reason to think that tho tarzet LiF-1, which was used 
for some of the earlier data is somewhat uneven. Since the beam 
without the foil interposed is concentrated bit is more diffuse with 
the foil interposed, the portion of spread in corresponding peaks due 
to tarset thiclmess may not cance! in the subtraction of the peak 
wer2ances. 

(b) As indicated in Section III.B, very fine, undetectable thick- 
ness variations nay cause a stra:slins efféct in addition to the 
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V, Summary of Results. 


The results on stopping power are summarized in Table XXXV, 
giving the values for stopping power in kev=cn* ng, and in 
Table XXXVI which gives the results for atomic stopping power ( O^) 
in ev-en^, Figure 11 gives in sraphical form the result for 
stopping power, while Figures 12, 13, and 14 give in graphical 
form the atomic stopping power results expressed as 9, jt? 

The results in stopping pover straggling, e:xyressed as 


= (n°) are given in Table XXXVII and Figure 15. 


57.4 1m 





Table XXXV. (Two Pages) 


SUMMARY OF DATA ON STOPPING POWER FOR PROTONS 


(kev-en* /mg) 





(a) Nickel | (b) Copper 
Proton Stopping prob, Proton Stonving Prob. 
Energy Power Frror, Energy Power Error, 
(xev) % (kev) % 
527 172.0 1.5 446 169,0 1.5 
704 148,7 2.0 532 153,0 1.5 
718 150.0 5.0 603 150.1 1.5 
739 142.8 2.0 713 139.9 1.29 
741 141.7 2.0 755 137.0 2.0 
755 145.3 1.5 812 151.6 1.5 
757 146.3 . 5 949 120,0 2.0 
915 18 30 996 114,7 xb 
935 125.7 2.0 1006 112.0 2.0 
941 133.7 2.0 1050 119.1 2.0 
949 128.9 2.0 
951 121.5 2.0 =н να 
977 129,0 5.0 
1000 127.4 2.0 
1007 121.9 2.0 (а) Nitrogen (#No) 
1046 170,4 2.0 
1047 122,0 2,3 Proton Stopping Prob. 
1057 120.7 2.0 Energy Power Error, 
(kev) % 
406 400 4.0 
c) Argon 409 401 4.0 
452 371 2.5 
Proton Stopping Prob, 501 362 4.0 
Energy Power Error, 502 351 4.0 
(kev) e 54] 328 5.0 
| 648 278 4.0 
421 224.5 3.9 648 296 4,5 
476 224 2.0 716 275 4.5 
519 215 3.5 717 264 4.5 
567 207 2.0 748 297 3.0 
662 176 4.0 914 219 4,5 
713 187 200 914 238 4.5 
778 176 2.5 942 ` 217 3.9 
929 154 4.0 973 207 4.5 
932 155 4.0 974 228 4.5 
971 154 5.0 1001 209 4.0 
989 148.5 4.0 


Table XXXV. (Concl.) 


(e) Neon (£) Krypton 

Proton Stopping Prob. Proton Stopping Prob. 

Energy Power Error, Energy Power Error, 
(kev) % (kev) 6. 
422 310 3.9 420 163 3.5 
424 290 3.5 516 149,5 3.5 
517 280 3.9 527 153.5 3.5 
519 269 5.0 665 130 4,0 
668 235.5 4.0 675 129 5.0 
730 217 4.0 733 124 4,0 
733 217 4.0 743 122 3.5 
929 186 4.0 929 107 4.0 
930 191 4.0 941 106.5 4.0 
988 175 4,0 989 102.5 4.0 
992 188 4.0 

(g) Xenon 

Proton Stopping Prob. 

Energy Power Error, 
(kev) % 
457 120.5 3.5 
441 121 5.0 
536 110.5 3.5 
938 115 3.5 
679 91 3.5 
679 95 3.5 
748 86,9 3.5 
750 92.7 2.9 
943 77.8 4.0 
969 77 3.5 


SUMMARY OF DATA ON ATOMIC 


STOPPING POWER FOR PROTONS (10715 ew-em”) 





(a) Nickel b) Copper 
Proton Atomic Prob. Proton Atomic Prob, 
Energy Stopping Error, Energy Stopping Error, 
(Kev) Power % (Kev) Power % 
527 16.78 1.5 446 17.86 1.5 
704 14,48 2,0 552 16.70 1.5 
718 14.62 5,0 603 15.87 1,5 
739 13.91 2.0 713 14.78 1.5 
741 13.80 2.0 755 14.48 2.0 
755 14.18 1.5 812 13.90 1.5 
757 14.27 1.5 949 12.68 2,0 
915 | 12,50 3,0 996 12.11 1.5 
935 12.85 2.0 1006 11.83 20 
941 13,02 20 1050 11.95 2.0 
949 12.51 2.0 
951 12.51 2.0 
977 12.57 3.0 
1000 12.42 2.0 
1007 11.88 2.0 
1046 11.74 2.5 (a) Nitrogen (3N5) 
1047 11.89 2.9 
1057 11.77 2.0 Proton Atomic Prob, 
Energy Stopping Error, 
(Kev) Power % 
c) Argon 406 9.30 4,0 
409 9.32 4,0 
Proton Atomic Prob. 452 8.62 2.5 
Energy Stopping Error, 501 8.42 4.0 
(Kev) Power % 502 8.16 4,0 
541 7,62 3.0 
421 14.87 3.9 648 6.47 4.0 
476 14.81 2.0 648 6.87 4.5 
519 14.28 3.5 716 6,40 4.5 
567 13.71 2.5 717 6.15 49 
662 11.67 4.0 748 9. 98 3.0 
713 12,39 Leo 914 5.09 4,5 
778 11,66 2.9 914 9.93 4.5 
929 10.19 40 942 9.05 38 
932 10.24 4.0 973 4.80 4,5 
971 10.20 3.0 974 5.50 4.5 
989 9,82 4.0 1001 4,8 4,0 
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e) Neon f ton 





Proton Atomic Prob, Proton Atonie Prob. 

Energy Stopping Error, Energy Stopping Error, 
(Kev) Power % (Kev) Power % 
422 10,37 3.9 420 22. 60 Jeo 
424 9.70 5.90 516 20.75 5.9 
517 9.38 9.0 527 21.30 3.5 
519 8.99 5.9 665 18.07 4.0 
668 7.88 4.0 675 17.95 3.5 
750 7.27 4.0 733 17.21 4.0 
755 7.26 ' 4,0 743 16.98 5.5 
929 6.22 4.0 929 14.81 4.0 
950 6.58 4,0 941 14.77 4.0 
988 5.86 4.0 989 14.22 4.0 
992 6.50 4.0 

( g) Xenon 

Proton Atomic Prob. 

Energy Stopping Error, 
Kev Power 
437 26.30 3.5 
44] 26.40 3.0 
536 24.10 3.9 
558 25.10 3.9 
679 19.85 5.9 
679 20.72 Jed 
748 18.95 3.5 
750 20.22 8.9 
943 16.99 4.0 
969 16,80 de 9 








Table XXXVII. 


SULMARY OF DATA ON STOFPING PO ER STRAGGLING OF 


PROTONS ( = (Ω2)) 


(kev?-cm* /ng) 
(a) Nickel 
Proton Energy Stragsling in Stovning Tower 

(kev) (Varience per Unit Thickness) 
527 88.8 

704 77.6 

739 67.6 

741 61,4 
757 53.9 

941 94.0 

1046 116.2 

(Ъ) Copper 
Proton Energy Straggling in Stopping Power 

(kev) (Variance per Unit Thickness} 
446 34.2 

749 62.8 

756 56.5 

945 64.9 

949 72.0 

1004 75.6 

1006 95.5 
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VI, Conclusion ~ Discussion of Results. 


A. Comparison of Present Results with Other Published Results. 


1. Stoppin” Power = ison with Theory 





In Figure ll, we have plotted mean curves obtained from our 
results on atomic stopping power. Superimposed on our resulting 
curves are curves obtained from theoretical predictions. The set 
of curves on the right side of the figure are obtained hy use of 
equation (I-B-69), as modified ly inclusion of the correction Cy 
in the mamer indicated by equation (I-B-63). The values of I y 
were obtained from Mano (1934), and are listed in Table II. The 
values of C, were obtained from the graphs of Walske (1952). As 
this theoretical computation camot be expected to be very accurate 
in the energy region of interest to us , (primarily due to the 
failure to include corrections for deficiency in stopping power of 
L and hizher shell electrons), it is not surprising that the theo- 
retical curves tend toward too high a value as the proton energy 
decreases. It is seen that a smooth curve of transition could easily 
be dram fran our experimental curves to approach the theoretical 
curves at higher energies, 

The comp:ted curves of Hirschfelder and Magee (1948), which 
traverse the energy region of interest to us, are plotted for the 
gases argon and xenon, Since their approach has included the stopping 
power corrections for all shells, one may expect their curves to 
agree better with experimental data than those based on the older 


formulations. In both cases the experimental curves lie below phe 
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theoretical curves by about 5 to 10%. It would be unreasonable to 
expect better agreement (except by accident), since the theorctical 
considerations on which Hirschfelder and Magee based their analyses 
were admittedly approximate. 

Livingston and Bethe's (1937) curve for proton atomic stopping 
power in air is likewise plotted, as a rough comparison with our 
results for nitrogen, Here, too, experiment seems to be about 10% 
below the predictions in this energy region, 

No effort is made to include the predictions by Bohr's (1948) 

formula for slow particles in heavy media (Equation I-B-24), since 
it is not expected to give quantitatively correct answers. Asa 


matter of fact, the formula gives predictions for 0 Ja /3 of 


2 
8.5 ev=cm at 500 kev proton energy and 6.1 even” at 100 kev. 


Comparison with Figure 14 shows that these values are consistently 
higher than experiment by an appreciable amount. The theory can 

be checked however for_correctness of dependence on Z and v. 

The plots of σ-/αν all seem to show thet for heavy media the 
value of this variable is almost independent of Z This conforms 
to the predicted variation of atomic stopping power with atomic 
number - at least for elements of atomic number >30, The in- 
verse dependence on v predicted by Bohr's theory can be checked 
by determining from Figure 11 whether the dependence on E is as 
the inverse square root of the latter variable. Since Figure 11 

is plotted on a log-log scale, the slope of the curve should deter- 
mine this exponent. It is seen that all the curves have a slope of 


approximately minus one-half, in accordance with the theory. The 
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heavier cases ore more nearly in accord, as is expected. The lighter 

sases show a steeper slope, and the metals a somewhat less steep 

slope, the anomalous results from the netals will be discussed 

in the concluding section, 

2. Stopping Power -ü Comparison with Other Experimental Information, 
The sources of other experinental information have already been 

indicated in Section I,C. Figures 12, 13, and 14 give data on 


er pal 3 


in such a way that the comparison of our results with that 
reported from other laboratories is readily seen. 

In regard to the metals, Firure 12 includes results, attributed 
to Kahn. Our curve for copper is about 4% below that of Kahn, 
which discrepancy is not alarming. Kahn has matched his curves to 
Warshaw's (1949) work at the 350-400 kev region, but our curve for 
copper can be matched with Warshaw's just as easily. Тһе data of 
Madsen (1953) on copper is much different from other experimental 
data and is not included. The earlier work of Madsen and Venkates-~ 
warlu (194£a) on beryllium is essentially in — vith the 
curve fiven by Kahn, so that a separate curve for this contribution 
is unnecessary im Figure 12, although credit is certainly due to 
these earlier workers, 

Our data for gases is compared in Fisure 13 with the data fron 
the California Institute of Technology up to 600 kev, reported by 
Dunbar, et al. (1952). The work of Weyl (1952) on the lighter gases, 
especially on argon, up to the maximum ener Ύ he utilized of 400 kev, 


agrees with the data from the California Institute of Techmolosy, so 


thet his results are not explicitly included. It cen be seen tht 


πι κ 


our data agree well with other results for the lighter gases. The 
data for krypton differ by a little more than the probable errar, 
but the results for xenon differ by much more than ihe proboble 


error, 





. Comparison of Results of Strazzling Measurements with Theory. 
There are no experimental measurements of stopping power 
straggling with which we can campare our results on nickel and 
copper. Since Madsen and Venkateswarlu's (1948b) measurements 
for beryllium checked so well with the theory of Bohr (1942), 
it is of special interest to check our results with theory. 

Equations (I-B-95) and (I-B-96) provide the relation, accord- 
ing to Bohr (1948), for theorctieal prediction of straggling, 
They are readily evaluated in the cases of the metals nickel 
and copper. We should note that according to (I-B-96) the number 
of effective electrons, 2 » becomes about 27 for 500 kev protons 
and 29 for 1000 kev protons. Since the latter figure is considerably 
nore than the total number of electrons in the copper or nickel atom, 
it is obvious that the estimate is too high. We have seen in sub- 
section 1, however, that Bohr's use of the Thomas-Fermi model gives 
too high an answer for stoppning porer, which leads us to the suspi 
cion that the formula obtaincd from the model (I-B-72) gives too 
large a valuc. 

Evaluation of the theoretical predictions according to Livingston 
and Bethe (1937) is a little more difficult. (Sce equation I-B-97). 


The values of I! are not known precisely, nor is IC. easy to evolu- 
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ate. A calculation was made indicating the contribution made by the 
term in which these constants appear, using 4/3 for each I n and 
the ionization energy I, in place of the average excitation ener- 


& Il. This gives a value to this expression of about 8.5 for 

500 kev protons and about 11.5 for 1000 kev protons, valid both for 
copper and nickel. Гог 42! we mst subtract from the atomic numbers 
of the elenent studii the oscillator strength of the K and L electrons, 
in accordance with the criterion of equation (I-D-71). Using the 
values computed by Honl (1933) for oscillator strength we get for 

4' values of 19,8 for nickel and 20,8 for copper. 

The theoretical curves and experimental data are plotted in 
Figure 15. It can be seen that our reais are too scattered to per- 
mit us to make ay positive statements. The following conclusions seem 
to be indicated, however: The theory of Livingston and Bethe is in- 
perfect in that it does not permit gradual changes in the value of 
4! , the effective number of electrons; that of Bohr takes account 
of this change with proton energy, but does not include the extra 
term sterming from the kinetic energy effect of theq@bital electrons. 
On the other hand, the equation of Bohr over-estimtes the number of 
effective electrons, and in this particular rezion this error nearly 
cancels the failure to include the kinetic energy effects. Tie ex- 
perimental results roughly follow the trends with energy as predicted 
by Bohr (1948) thagh they seem to indicate that in absolute value, 


theory may give too high a prediction in this region. 





By Concluding Comments. 


This dissertation should not be concluded without sore comment 
on the general signifiance of the experimental results on stopping 
pover in this energy region, There are a fev discrepancies between 
the results reparted by the various experimenters, including the 
writer, but nevertheless there is sufficient agreement to detect the 
pattern of the stoppins power effect at these moderately low proton 
energies. Figures 12, 13, and 14 - especially the latter - bring 
out a very interesting fact. Above 700 kev there is little differ- 
ence between the stoppinz power of gases and the solids studied ex- 
cept the predicted variation with atomic number. Below 700 kev 
however the results indiatte a consistency anong the gases or among 
the metals, tit a definite variation between the two groups of ele- 
ment s. 

Theory has not developed to the eztent of permitting a definite 
explanation of this effect. In a qualitative way, however, we may 
speculate upon the possibility that those special effects in con- 
densed media (condutting and non-conducting) discussed in Section 
I.B.6 may provide the explanation desired. Whereas the polariza- 
tin effects in cmdensed media of hich atomic number are considered 
to have a negligible effect in view of the sall? proportion of 
loosely bound electrons contributing to polarization, this presunp- 
tion may no longer be valid at low proton energies at which only 
the outer electrons contribute to the stopping power phenomenon, 


Since theoreticians have admitted that deficiency in stopping nay 
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exist for both cnnducting and non-conducting types of condensed nedia, 
we cannot state whether our resulting difference between the sases 
and metals in stopping of bw enerzy protons is due to the special 
conducting properties of the metals, or simply to their condensed 
state. It is proper and desirable, of course, to make no final judg- 
ment on sich a matter until experiments are forthcoming on the 
stopping power of non-conducting elements in the solid state far pro= 
tons in this enerzy region. 

One eunent nay be made on the stopping power of neon, which 
according to Fimre 14 appears to have an anomalously low value of 
stoppin; power, This is not too surprising however upon sone con- 
sideration of theory. The theoretical prediction of smooth vari- 
ation with atomic number is a result of Bohr's use of the Thomas- 
Fermi statistical model of the atom, which is well-known to be 
least valid for the lisnt elements. Since neon contains only closed 
shells, its outer eleebrons lie deeper in the Coulombian potentiel 
well than those of the atoms adjacent in the atomic table. Its 
electrons are therefore likely to be less effective thm those of 
the neighboring atoms, even aside from the general trend with atomic 
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Theoretícal Triode Characteristic Curves 


Figure App-2 


Figure Appel shows the basic circuit of the voltage regulator 
for the ramet circuit. Figure App=2 gives the characteristic curves 
for an ideal triode. The tube used, GAS7, has characteristic curves 
quite close to the ideal case = at least in region in which we expect 
to operate. The characteristics may be theoretically expressed 


T P 





by the following formula: 


b. . — ; (App-1-1) 


where AC is the tube amplification factor, R. is the plabe resis- 
tance of the tube, and the other symbols are obvious from the 
above figures. 


Fron Figure Арр-1 we can see that 


E wIRA S , (App-1-2) 
and 

E = F - (I+I)(R +R) . (Anp-1-3) 
Also we note that 

wir ΠΤ} SUNL (App- I-A) 


The above four independent equations can be reduced to a single 
equation by the elimination of the three variables Е, ; E», and T, š 
By doing so we may get an equation for I , which can be reduced to 


the following forms 


τ - A _ y (App-I-5) 
R t R. +Z) +— 


R, +R, 
R, = AB, ° (App-1-6) 


From these last two equations we note that I is independent 


with 
r 


of ES when г = O , or when R £ Ry/M . Sinee this ls 
factor is the reciprocal of the transconductance of tho tube, the 
criterion that output current (and output voltage) be constant is that 
conductance of the plate resistor equal the transconductance of the 


tube = whence the name of the regulator type: "transcond::ctancoe", 





Unfortunately no actual tube characteristic curves correspond 
to the ideal curves, so that different operating points on the 
actual chart have different values of transconductance, The proper 
value for Ry is taken then as the reciprocal of triode transcon~ 
ductance at some mean operating point located in the most linear 
portion of the actual set of curves. 

Once R. is selected then r nay vary from zero only if the 
characteristic point at which we are operating is one with a dif- 
ferent transconductance value than the one used to select Rm. I 
the actual characteristic curves are not too different from the 
ideal (this is one reason for selection of the 6AS7 tube — the 
others are its high current capacity and its low transconductance, 
permitting a reasonably low value for R.) ; then r will not deviate 
markedly from zero. We see thus from equation (Apn-I-5) that I 
is affected very little by variations in r of the order of mag- 
nitude expected. 

We note also that if r is approximately zero, I is not only 
practically independent of E, but also of variations in R, . Ths 
resistance is not needed theoretically but has a very good use in 
actual practice. Assuming that r is small enough to be considered 
zero, we may use equations (App-1-2), (App-1-3), and (App-I-5) to 


determino that: Р 


А 1 d (Rp Rj FRS)(A AV.) 
l7 (m + E) a R + R (+4) 


(App-1-7) 
We see that In does depend simificantly on πο « This enables 


us to vary 5 appreciably without markedly changing I or E, ο 
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The value of this will be seen in a monent, 

Ihe systei must be flexible enough so that any desired value of 
regulated current may be obtained from about 0.4 amperes to about 1.2 
amperes, (During the process of cycling the magnet, a maximum 
current of not less than 3 amperes is desirable, but precise regu= 
lation is not required in this process, and the desim is made so that 
the resulator may be removed from the system by a single switching 
operation.) This means that E, rust be allowed to vary at will 
without affecting the regulation effect, It happens that lines of 
constant transconductance can be drawn on the charactoristic chart 
for the tube selected, which lines proceed in general from the 
lower left side of the chart to the upper right side. If the 
system is designed so that the operating point, even if it changes, 
stays on or near the line selected to match Ry ə then the regula- 
tion property renains for all operating currents, In case the sys- 
tem characteristics cannot be nade to follow the proper transconduc- 
tance line (which we nay call the operating line) exactly, changes 
in the value of R, can be nade to bring one exactly to said opera- 
ting line, In practice it is quite easy to determine the proper value 
of Ro to keep the system operating at the proper transconductance, 
With r slightly different from zero, any chanze in R, will make 


2 
a just perceptible change in current I. R, is simply changed until 
the position is reached for which slight changes make no appreciable 
chances in mamet current, 
Gross chanzes in magnet cirrent are effected by chances in E, > 


which is caused in turn by changing the setting on the field current 


SIF - 





variac, as explained in Section III.C. In order to prevent the 
regulator from opposing this desired current change, it is necessary 
to change simultaneously some other factor in the equation (App-I-5) 
on which the current depends. Obviously Vo is the only one which 
it is practicable to change. We provide this bias by means of a 
bank of batteries across 2 potentiometer, the middle connection of 
which, being variable, can provide a variable value of Ve e Taig 
potentiometer is then "ranged" with the variac, so that they are 
changed together. By proper selection of the potentiometer resistance 
and battery voltage, the value of prid bias will change along with 
the plate voltage (E,) so that the operatin$g point moves approxi- 
mately along the line of proper transconductance on the characteris- 
tic chart, 

The complete diagram of the magnet control panel, including 
all the control and regulator elements, is given in Figures App=3(a) 
and App~3(b), separated into the upper and lower sections, This 
separation can easily be carried out in practice if necessary, 


permitted by the use of multiply-canged connectors betuocen the two 


sections, 
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